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DISPOSITIVO E METODO DI FILTRAGGIO PER LA RIDUZIONE DEL 
RUMORE IN SEGNALI ELETTRICI, IN PART ICOLARE SEGNALI 
ACUSTICI E IMMAGINI 

5 La presente invenzione riguarda un dispositivo ed 

un metodo di filtraggio per la riduzione del rumore in 
segnali elettrici, in particolare segnali acustici (vo- 
cali) e immagini . 

Come e noto, per rimuovere component i indesiderate 

10 da segnali elettrici sono utilizzati filtri operant! 
secondo diverse tecniche lineari e non lineari. In par- 
ticolare, le componenti - indesiderate possono essere un 
qualsiasi tipo di rumore (bianco, flicker, .etc.) o al- 
tri tipi di segnali acustici o visivi sovrapposti. 

15 I filtri ■ lineari -sono attualmente lo strumento piu 

utilizzato per il filtraggio del rumore. I filtri FIR 
(dall ' inglese Finite Impulse Response) eliminano tutte 
le armoniche di un segnale che hanno una frequenza mag- 
giore della frequenza di taglio del filtro e migliorano 

20 il rapporto segnale/rumore (SNR) . Un'altra tecnica di 
filtraggio lineare e basata sulla trasformata veloce di 
Fourier (FFT-Fast Fourier Transform) , in cui il segnale 
viene trasf ormato nel dominio delle frequenze, le armo- 
niche indesiderate sono rimosse, e quindi viene calco- 

25 lata la trasformata di Fourier inversa. 
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Per quanto riguarda le tecniche non linear!, fil- 
tri median! sono semplici da progettare e consentono di 
realizzare circuit! hardware semplici. I filtri median! 
si basano sul confronto dei singoli campion! di segnale 
in un intervallo di tempo con la media di tutti i cam- 
pion! nello stesso intervallo di tempo. In base a tale 
confronto, i singoli campion! sono attenuati selettiva- 
iaente. 

Tutti questi metodi hanno lo svantaggio che, du- 
rante la rimozione del rumore, vengono rimosse anche 
alcune componenti del segnale originario. 

Inoltre, nessuna delle attuali tecniche, sia li- 
neare che non lineare, compreso il filtraggio mediano, 
e in grado di preservare front! ripidi del segnale. Se 
si usa un filtro a media mobile ("moving average"), 
1'ampiezza della finestra deve essere molto piccola se 
si desidera preservare front! ripidi. Tuttavia, se la 
dimensione della finestra diventa piccola, non si ha 
una riduzione signif icativa dell ' energia del rumore. Se 
si usano filtri lineari, tutte le frequenze al di sopra 
della frequenza di taglio vengono eliminate, con conse- 
guente forte distorsione del segnale. 

Scopo dell'invenzione e quindi mettere a disposi- 
zione un metodo ed un dispositive di filtraggio senza 
causare un sensibile deterioramento del segnale e pre- 



129-11-2000] 



servando i fronti del segnale stesso. 

Secondo la presente invenzione vengono realizzati 
un dispositivo ed un metodo di filtraggio di segnali 
elettrici, come definiti nelle rivendicazioni 1 e, ri- 
5 spettivamente, 24. 

II filtro e il metodo descritti sono particolar- 
mente utili nel caso di segnali aventi fronti ripidi, 
per i quali si desidera preservare i fronti dei segna- 
li. Inoltre, e possibile filtrare segnali affetti da 

10 rumore bianco e non, quale ad esempio rumore flicker. 
Grazie al presente metodo e inoltre possibile eliminare 
da un segnale altri segnali sovrapposti ad esso con ca- 
ratterizzazione spettrale svariata. 

II dispositivo e il metodo descritti sono basati 

15 ■ su una rete neuro-fuzzy. Essi sono. implementati con" una 
tecnica di filtraggio a media mobile in cui i pesi per 
la ricostruzione finale del segnale sono calcolati in 
una rete neuro-fuzzy in base a specifiche regole fuzzy. 
Cio consente una migliore riduzione del rumore. Le re- 

20 gole fuzzy operano su diverse variabili, chiamate ca- 
ratteristiche ("features") di segnale. Qui sono de- 
scritte tre caratteristiche di segnale e sei regole 
fuzzy. 

Il filtro proposto e adatto a segnali visivi o se- 
25 gnali acustici, anche con variazioni improvvise. Per 
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creare le regole possono essere utilizzati vari tipi di 
funzioni o caratteristiche di segnale. Con il metodo 
descritto, le caratteristiche di segnale sono correlate 
alia posizione del campione nella finestra di campioni 
5 considerata, alia differenza fra un campione e il cam- 
pione al centro della finestra e alia differenza tra un 
campione e la media di campioni nella finestra. Queste 
caratteristiche di segnale possono influenzare notevol- 
mente i valori dei pesi per la ricostruzione del segna- 
10 le; inoltre esse possono essere calcolate in maniera 
relativamente semplice. 

II metodo e il filtro secondo l'invenzione com- 
prendono inoltre un banco di filtri neuro-fuzzy. In tal 
modo, il segnale pud essere separato in diverse sotto- 
15 bande tramite la teoria delle ondine ("wavelets"); cia- 
scun segnale di sotto-banda puo essere filtrato da una 
rete neuro-fuzzy e quindi le varie sotto-bande possono 
essere ricostruite dal banco di filtri di sintesi. Come 
e noto dalla teoria "wavelet", nella prima sotto-banda 
20 le caratteristiche di segnale hanno bassa frequenza, 
mentre nell' ultima sotto-banda le caratteristiche del 
segnale hanno la massima frequenza. Se deve essere can- 
cellato del rumore non bianco (rumore colorato) , questo 
viene approssimato da rumore bianco in ciascuna singola 
25 sotto-banda. Dato che una rete neuro-fuzzy lavora bene 
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sul rumore bianco, questa soluzione porta ad una sensi- 
bile riduzione del rumore. 

La rete viene addestrata fornendo alcune configu- 
razioni di segnali di ingresso di uscita (la configura- 
zione del segnale di uscita che si vuole ottenere per 
effetto dell ' evoluzione della rete e chiamata configu- 
razione desiderata o "target"). L'algoritmo di adde- 
stramento e basato su uno dei metodi di apprendiraento 
noti, quali il "gradient descent", un algoritmo geneti- 
co, il metodo di "annealing" simulato, la ricerca ,ca- 
suale o qualunque altro metodo per 1 ' ottimizzazione . di 
f unzioni . 

Per la comprensione della presente invenzione ne 
vengono ora descritte forme di realizzazione preferite, 
a puro titolo di esempio non limitativo, con riferimen- 
to- ai disegni -allegati, nei quali: 

- la figura 1 mostra uno schema a passi mostrante 
1 T architettura generale del filtro secondo 1' invenzio- 
ne ; 

- la figura 2 rappresenta la topologia di un passo 
di figura 1 relativo ad una rete neuro-fuzzy; 

- le figure 3a e 3b mostrano uno schema di flusso 
relativo alle operazioni eseguite dal filtro di figura 
1; 



- la figura 4 mostra uno schema a passi di un di- 




spositivo di filtraggio utilizzante il filtro di figura 
1; 

- la figura 5 mostra le funzioni di trasf erimento 
di filtri digitali appartenenti al dispositivo di fil- 

5 traggio di figura 4; 

- la figura 6 mostra il diagramma a blocchi di al- 
cuni filtri digitali appartenenti al dispositivo di 
filtraggio di figura 4; 

- la figura 7 mostra uno schema a passi di un al- 
io tro dispositivo di filtraggio utilizzante il filtro di 

figura 1; 

- la figura 8 mostra uno schema di flusso relativo 
al funzionamento del dispositivo di filtro di figura 7; 

- la figura 9 mostra uno schema a passi relativo 
15 all ' implementazione hardware del filtro di figura 1; e 

- le figure 10a, 10b, 10c mostrano rispettivamente 
un segnale vocale privo di rumore, lo stesso segnale 
sovrapposto a rumore bianco e lo stesso segnale dopo il 
filtraggio secondo 1 ' invenzione . 

20 La figura 1 mostra un filtro 1 comprendente 

un'unita di calcolo caratteristiche segnale 2, una rete 
neuro-fuzzy 3, un'unita di ricostruzione 4 e una unita 
di addestramento 5. 

L'unita di calcolo caratteristiche segnale 2 rice- 

25 ve in ingresso un segnale In formato da una pluralita 
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di campioni di ingresso e(i), memorizza, ad ogni ciclo, 
(2N+1) campioni di ingresso e(i) (che rappresentano una 
finestra di lavoro per il filtro 1) in un buffer inter- 
no, calcola le caratteristiche di segnale XI (i), X2 (i) 
5 e X3(i) per ogni campione di ingresso e(i) sulla base 
di tutti i campioni di ingresso appartenenti alia fine- 
stra di iavoro (come descritto in dettaglio in seguito) 
e fornisce le caratteristiche di segnale XI (i), X2(i) e 
X3(i) cosi calcolate alia rete neuro-f uzzy . 3 . 

10 La rete neuro-fuzzy 3 elabora le caratteristiche 

di segnale XI (i) , X2(i) e X3(i) e genera in uscita un 
peso di ricostruzione oL3(i) per ogni campione di in- 
gresso e(i). A tale scopo, a partire dalle caratteri- 
stiche di segnale XI (i) , X2(i) e X3(i) e per ogni cam- 

15- pione di ingresso e(i).,. la rete neuro-fuzzy 3 esegue 
dapprima un ' operazione di fuzzificazione," quindi appli- 
ca regole fuzzy prefissate e infine esegue un'operazio- 
ne di def uzzif icazione . II peso di ricostruzione oL3 (i) 
cosi ottenuto e quindi la somma pesata di tutti i cam- 

20 pioni di ingresso e(i) nella stessa finestra di lavoro, 
come spiegato in dettaglio in seguito con riferimento 
alle figure 2 e 3. 

L'unita di ricostruzione 4 riceve i pesi di rico- 
struzione oL3(i) e i campioni di ingresso e(i) e, dopo 

25 aver accumulato un sufficiente numero di campioni di 



ingresso e(i) e di relativi pesi di ricostruzione 
oL3 (i) , genera un campione di uscita u(i) la cui suc- 
cessione forma un segnale di uscita Out. 

L'unita di addestramento 5 e operante solo ini- 
5 zialmente, in modo da addestrare la rete neuro-fuzzy 3 
e modificare i pesi della rete alio scopo di ottenere 
un comportamento ottimale del filtro 1, come sotto de- 
scritto in dettaglio. 

Le caratteristiche di segnale XI (i) , X2(i) e 
10 X3(i), calcolate nell'unita di calcolo caratteristiche 
segnale 2, sono correlate, rispettivamente, alia di- 
stanza fra ciascun campione e il campione centrale del- 
la finestra considerata, alia differenza tra un campio- 
ne e il campione al centro della finestra, e alia dif- 
15 ferenza tra un campione e la media di campioni nella 
finestra stessa, e sono normalizzate in modo da ottene- 
re valori compresi fra 0 e 1. 

In dettaglio, data una finestra di (2N+1) campioni 
di ingresso e(i), con i = 0, 2N, le caratteristiche 

20 di segnale XI (i) , X2(i) e X3(i) per ogni campione di 
ingresso e(i) sono definite come: 

max.(aijj ) 
\e(f)-av\ 

AT3(i)= 1 1 : (3) 

msoi(diff _av) 
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in cui N e la posizione di un campione centrale 
e (N) nella finestra di lavoro; 

max (diff) = max { e (k) -e (N) } con k = 0, 2N, cioe 

il massimo delle differenze fra tutti i campioni di in- 
gresso e(k) e il campione centrale e (N) ; 

av e il valore medio dei campioni di ingresso 

e(i); 

max (diff_av) = max {e (k) -av} con k = 0,..., 2N, cioe 
il massimo delle differenze fra tutti i campioni di in- 
gresso e(k) e il valore medio av. 

La rete neuro-fuzzy 3 e una rete fuzzy a tre stra- 
ti la cui rappresentazione funzionale e mostrata in fi- 
gura 2, nella guale, per semplicita, 1 1 indice i tra pa- 
rentesi relativo alio specif ico campione all ' interno 
della rispettiva finestra. di lavoro non e indicato . . 'Cio. 
nonostante, come sopra indicato e come risultera piu 
chiaro dallo schema di flusso delle figure 3a, 3b, 
1 ' elaborazione neuro-fuzzy rappresentata dalla figura 2 
viene ripetuta per ogni campione di ingresso e(i). 

In dettaglio, a partire dalle tre caratteristiche 
di segnale XI, X2 e X3 (o genericamente da 1 caratteri- 
stiche di segnale XI) e date k funzioni di appartenenza 
di tipo gaussiano per ogni caratteristica di segnale 
(descritte dal valore medio VI m (l,k) e dalla varianza 
W v (l,k)), viene eseguita una operazione di f uzzif icazio- 
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ne nella quale viene valutato il livello di appartenen- 
za delle caratteristiche di segnale XI, X2 e X3 rispet- 
to a ciascuna funzione di appartenenza (qui due per 
ogni caratteristica di segnale, per cui k=2; in totale, 
si hanno M = 1 x k = 6 funzioni di appartenenza) . 

In figura 2, tale operazione e rappresentata me- 
diante sei neuroni di priiao strato ("first layer neu- 
rons") 6 che, a partire dalle tre caratteristiche di 
segnale XI, X2 e X3 (genericamente indicate come XI) e 
usando come pesi il valore medio W m (l,k) e la varianza 
W v (l,k) delle funzioni di appartenenza, forniscono cia- 
scuno un'uscita di primo strato oLl(l,k) (in seguito 



indicata anche come oLl (m) ) calcolata come segue 



oIAQ, k) = oL\{m) - exp 



' r Xl-W m {l,k) \ 
I W v {l,k) J J 



(4) 



20 



In seguito, viene effettuata una operazione AND 
fuzzy utilizzando la norma del minimo, in modo da otte- 
nere N uscite di secondo strato oL2 (n) . N ad esempio e 
pari a 6. Come e noto, una operazione AND fuzzy utiliz- 
zando la norma del minimo si basa su regole del tipo: 
if X™ is and X? is A? and X? is ^ then 0 is 



if X™ is A[ n) and X? is A? and X^ is 4 n> then 0 is B\ 

in cui A?\ 4°, - B™, ecc. sono termini lin- 

guistic! quali "high" "low" e il valore dell'uscita O 
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per ogni regola e dato dal minimo dei livelli di appar- 
tenenza moltiplicato per un peso. 

In pratica, con la rete neuro-fuzzy di figura 2, 
ogni uscita di secondo strato oL2 (n) e pari al minimo 
5 fra i prodotti tra le M uscite oLl (m) dei neuroni di 
primo strato 6 ed un rispettivo peso di secondo strato 
W FA (m,n) . 

In figura 2, tale operazione e rappresentata me- 
diante N neuroni di secondo strato 7 che implementano 
10 l'equazione: 

oL2(n) = rmn{tV FA (m,ri)-oL\(m)} (5) . 

n 

Infine, il terzo strato corrisponde ad un' opera- 
zione di def uzzi f icazione e fornisce in uscita un peso 
di ricostruzione oL3 di tipo discreto, utilizzando N 
15 pesi di terzo strato W DF (n). II. metodo di defuzzifica- 
zione e quello del centro di gravita ed e rappresentato 
in figura 2 mediante un neurone di terzo strato 8 for- 
nente il peso di ricostruzione oL3 secondo l'equazione: 

•YW DF (ri)oL2(ti) 
oZ3=^-^ (6) 

n=\ 

20 L'unita di ricostruzione 4 attende quindi un suf- 

ficiente numero di- campioni e(i) e di relativi pesi di 
ricostruzione oL3 (1) (almeno 2N+1, pari alia larghezza 
di una finestra di lavoro) e calcola un campione di 
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uscita u(i) come somma pesata dei campion! di ingresso 
e(i-j), con j = 0 ... 2N, usando i pesi di ricostruzione 
oL3(±-j), secondo l'equazione: 
^oL3(i-j)-e(i-j) 

m-^t-i* (7) 

5 L'unita di addestramento 5 opera soltanto in una 

fase iniziale di apprendimento, durante la quale un se- 
gnale di ingresso In avente una conf igurazione nota 
viene fornito al filtro 1, il segnale di uscita Out ot- 
tenuto viene confrontato con un segnale desiderate (o 
10 di "target") Tg e la distanza fra il segnale ottenuto 
Out e il segnale desiderato Tg viene valutata sulla ba- 
se di una funzione di conformita ("fitness"). Tale fun- 
zione di conformita puo essere ad esempio il rapporto 
segnale/rumore quadratico avente espressione: 

in cui T rappresenta il numero totale di campioni 

di ingresso e(i). 

Utilizzando la funzione di conformita, e applican- 
do un metodo per 1 < ottimizzazione di funzioni, quale il 
20 "gradient descent", un algoritmo genetico, il metodo di 
"annealing" simulate, la ricerca casuale, vengono modi- 
ficati i pesi di primo strato (valore medio W n (l,k) e 
varianza W v (l,k) delle funzioni di appartenenza gaussia- 
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ne) e i pesi di terzo strato W DF (n) e viene eseguita una 
nuova valutazione di conformita. I pesi di secondo 
strato W FA (m,n) sono invece inizializzati in modo casua- 
le e non vengono modificati. II processo di apprendi- 
5 mento viene iterato fino a raggiungere un valore pre- 
fissato della funzione di conformita o fino a quando e 
trascorso un numero prefissato di generazioni dell'al- 
goritmo genetico o di passi dell 1 algoritmo di ottimiz- 
zazione selezionato. 

10 In questo modo la rete neuro-fuzzy 3 implementa, un 

algoritmo adattativo e supera le limitazioni delle reti 
neurali o dei sistemi fuzzy considerati separatamente . 
Infatti i sistemi fuzzy non hanno capacita di apprendi- 
mento e, se la scelta delle regole fuzzy non e acciira- 

15 ta, 1 ' algoritmo fuzzy non si comporta in modo soddis'fa- 
cente . Viceversa, utilizzando la rete neuro-fuzzy 3 ed 
eseguendo una fase di apprendimento preventiva, e pos- 
sibile approssimare un segnale avente una struttura ma- 
tematica complessa, ad esempio la voce in un ambiente 

20 rumoroso, senza conoscenza preventiva delle leggi mate- 
matiche che regolano tale sistema. 

II f unzionamento del filtro 1 di figura 1 e de- 
scritto in seguito in dettaglio con riferimento alle 
figure 3a, 3b. 

25 Inizialmente, il filtro 1 viene inizializzato, 



passo 10, cosi come un contatore di finestre p, passo 
11. Quindi il contatore di finestre p viene incrementa- 
to di (2N+1), pari al numero di campioni di ingresso 
e(i) in una finestra di lavoro, passo 12 e vengono ca- 
5 ricati i campioni di ingresso di una finestra di lavo- 
ro, compresi fra e (p) ed e (p+2N) , passo 13. Un contato- 
re di somme campioni h viene inizializzato a 0 ed una 
variabile di somma campioni S(0) viene inizializzata 
con il valore del campione di ingresso e (p) , passo 14. 
10 Quindi, il contatore di somme campioni h viene incre- 
mental di una unita, passo 15, e la variabile di somma 
campioni S (h) viene incrementata con il valore di un 
campione e (p+h) , passo 20. Dopo la somma di (2N+1) cam- 
pioni di ingresso, uscita SI dal passo 21, la variabile 
15 di somma campioni S (h) e pari alia somma di tutti i 
campioni di ingresso e(i) e puo essere direttamente 
utilizzata per il calcolo delle caratteristiche di se- 
gnale . 

In seguito, utilizzando le equazioni (1) , (2) e 
20 (3) , e il valore della variabile di somma campioni 
S(h), vengono calcolate le caratteristiche di segnale 
XI, X2, X3 per ogni campione e (p) -e (p+2N) , appartenente 
alia finestra di lavoro considerata, passo 22. Un con- 
tatore di campioni i, relativo ai campioni di ingresso 
25 all'interno della finestra di lavoro considerata, viene 
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quindi azzerato, passo 23, e successivamente incremen- 
tato di una unita, passo 24; un contatore di caratteri- 
stiche 1 relativo al numero di caratteristiche di se- 
gnale utilizzate (qui 3) viene azzerato, passo 25 e 
5 successivamente incrementato, passo 26; ed un contatore 
di funzioni k relativo al numero di funzioni di appar- 
tenenza utilizzate per ciascuna caratteristi.ca di se- 
gnale (qui 2) viene azzerato, passo 30, e successiva- 
mente incrementato, passo 31. Quindi vengono caricati i 
10 pesi di primo strato (il valore medio W m (l,k) e la va- 
rianza IV V (1 ,k) ) , passo 32, e vengono calcolate le usci.te 
di primo strato oLl(l,k), passo 33, in base all'equa- 
zione (4) . 

II caricamento dei pesi di primo strato W m (l r k), 
15 W v (l,k) e il calcolo delle uscite di primo strato 
o-Ll(l,k) vengono ripetuti per tutte le funzioni di ap- 
partenenza (uscita SI dal passo 34) e per tutte le ca- 
ratteristiche di segnale (uscita SI dal passo 35) . 

In seguito, un contatore di neuroni di secondo 
20 strato n viene azzerato, passo 39, e incrementato di 
una unita, passo 40; un contatore di algoritmo m viene 
azzerato, passo 41, e incrementato, passo 42; e vengono 
caricati i pesi di secondo strato W FA (m r n) relativi al- 
l'n-simo neurone di secondo strato 7, passi 43, 44. 
25 Quindi viene calcolata l'uscita di secondo strato 
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oL2(n) relativa all'n-simo neurone di secondo strato 7 
utilizzando l'equazione (5), passo 45, e la procedura 
dei passi 40-45 viene ripetuta per tutti gli n neuroni 
di secondo strato 7. 
5 Al termine (uscita SI dal passo 50), il contatore 

di neuroni di secondo strato n viene nuovamente azzera- 
to, passo 51, e increment ato di una unita, passo 52; 
vengono caricati in successione tutti i pesi di terzo 
strato W DF (n), passi 53, 54; al termine del caricamento 
10 (uscita SI dal passo 54), viene calcolato il peso di 
ricostruzione oL3(i) secondo l'equazione (6), passo 55. 

I pesi di ricostruzione oL3(i) cosi calcolati e i 
rispettivi campioni di ingresso e(i) vengono memorizza- 
ti, passo 60, e la procedura descritta dei passi 24-60 
15 viene ripetuta per successivi campioni di ingresso 
e(i), fino a raggiungere (2N+1) campioni di ingresso 
e(i) (uscita SI dal passo 61a). 

In seguito viene azzerato un contatore di soma j, 
passo 62, e il campione di ingresso e(i) viene sommato 
20 ai precedent! 2N campioni di ingresso e(i-j); la somma 
e memorizzata in una variabile E, passo 63. 

Dopo aver eseguito la somma di (2N+1) campioni di 
ingresso, viene calcolato un campione di uscita u(i), 
secondo l'equazione (7), passo 65. 
25 L'intero ciclo descritto dai passi 12-65 viene ri- 
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petuto fintanto che sono presenti campioni di ingresso 
e(i); al termine (uscita SI dal passo 70) termina la 
procedura di elaborazione dei campioni. 

Secondo un altro aspetto dell ' invenzione , il fil- 
5 traggio e basato su una analisi a multirisoluzione ot- 
tenuta mediante un banco di filtri in quadratura. La 
teoria delle ondine ("wavelets") fornisce la base teo- 
rica della multirisoluzione. 

Come e noto, una analisi a multirisoluzione defi- 

10 nisce un gruppo di sottospazi annidati di uno spazio 
funzionale sommabile in raodo quadratico ("square sum- 
mable function space", cioe lo spazio delle funzioni "di 
energia finita) , ampiamente noto nella fisica e nel- 
l'ingegneria elettrica. In base a tale analisi, una 

15 proiezione di una funzione in uno di questi 'sotto spa*zi 
isola i dettagli "piu grossolani" della funzione, men- 
tre la proiezione della funzione nel complemento orto- 
normale di questo sottospazio isola i dettagli "piu fi- 
ni" della funzione. La procedura puo essere iterata in 

20 modo da ottenere una piramide. Dalla teoria "wavelets" 
e noto che la stessa operazione puo essere eseguita 
tramite un banco di filtri FIR, in cui ciascun filtro 
FIR e seguito da un riduttore di campioni ("subsam- 
pler") . II segnale, cosi suddiviso in sottobande, puo 

25 essere ricostruito tramite un banco di filtri comple- 
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mentari, dotati ciascuno di un incrementatore di cam- 
pion! ( "upsampler" ) . 

Una prima forma di implement azione di tale solu- 
zione e mostrata in figura 4, ed e relativa ad un fil- 
tro multirisoluzione 80 generante due sottospazi e uti- 
lizzante una prima coppia di filtri FIR QMF (dall'in- 
glese Finite Impulse Response Quadrature Mirror Filter) 
per l'analisi ed una seconda coppia di filtri FIR QMF 

per la sintesi. 

In dettaglio, il filtro multirisoluzione 80 rice-' 
ve, su un ingresso 81, un segnale di ingresso s,(t). II 
segnale di ingresso s,(t) e fornito a due filtri FIR di 
ingresso H 0 , Hi che eseguono una operazione di convolu- 
zione in modo da fornire in uscita un rispettivo segna- 
15 le filtrato e 0 , e, pari a: 

('•)•<('->■) (9) 

in cui M e 1 ' ordine dei filtri H 0 , Hi; e 0 (t) , e,(0 e 
il t-esimo campione della rispettiva sequenza di usci- 
ta; s x (t) e il t-esimo campione della sequenza di in- 
gresso; ho(r), hi(r) e il tappo ("tap") t-esimo del fil- 
tro FIR di ingresso H 0 , Hi, in cui 
hi(r) = (-D r_1 h 0 (2I-r+l) 
con I numero intero . 
25 i filtri FIR di ingresso H 0 , H a hanno funzioni di 
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trasf erimento mostrate in figura 5 e definiscono quindi 
rispettivamente un filtro passa-basso ed un filtro pas- 
sa-alto ed hanno struttura mostrata in' figura 6, in cui 
il generico tappo h(r) corrisponde a h 0 (r) o a hi (r) , a 
seconda che si tratti del filtro FIR di ingresso H 0 o 
del filtro FIR di ingresso H-_ . 

Le uscite dei filtri FIR di ingresso H 0 , Hi sono 
collegate ciascuna ad una rispettiva unita di riduzione 
campioni { "subsampler " ) 84, 85 che scarta i campioni 
dispari dal segnale di ingresso e 0 (t) , e x (f) e mantiene 
solo i campioni pari, generando un rispettivo segnale 
e 0 (t), ei(t) . Le uscite delle unita di riduzione campio- 
ni 84, 85 sono collegate ciascuna ad un rispettivo fil- 
tro neuro-fuzzy 86, 87. I filtri neuro-fuzzy 86,. 87 
hanno • entrambi l'a struttura mostrata in figura 1 . 

I segnali di uscita u 0 (t), Ui(t) dei filtri neuro- 
fuzzy 86, 87 sono forniti ciascuno ad una rispettiva 
unita di incremento campioni ( "upsampler " ) 88, 89 che 
genera una rispettiva sequenza di uscita u 0u (t), Ui u (t) 
inserendo un campione zero fra ciascuna coppia di cam- 
pioni del rispettivo segnale di uscita u 0 (t), Ui(t) dei 
filtri neuro-fuzzy 86, 87. Le uscite delle unita di in- 
cremento campioni 88, 89 sono collegate ciascuna ad un 
rispettivo filtro FIR di uscita G 0 , G lf aventi anch'essi 
una rispettiva funzione di trasf erimento, date dalla 
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(9) e dalla (10) . 

Infine, i segnali di uscita dei filtri FIR di 
uscita Go, Gi sono sommati fra loro, campione per campio- 
ne, da un sommatore 92. 
5 La figura 7 mostra un filtro multirisoluzione 95 

utilizzante otto sottospazi . In dettaglio, il segnale 
di ingresso Si(t) viene inizialmente fornito a due primi 
filtri FIR di sintesi H 00 i e H 102 , rispettivamente di tipo 
passa-basso e passa-alto, e quindi viene ridotto di 
10 campioni ( "downsampled" ) in due prime unit A di riduzio- 
ne campioni 96, 97, analogamente a quanto descritto per 
le unita 85, 86 di figura 4. Le sequenze di campioni 
cosi ottenute sono fornite ciascuna a due filtri di 
sintesi (e quindi complessivamente a quattro secondi 
15 filtri FIR di sintesi Hon, Hu 2 , H013 e H n O • Le uscite 
dei secondi filtri FIR di sintesi Hon, Hm, H 0i 3 e H 114 
sono quindi fornite a quattro seconde unita di ridu- 
zione campioni 100-103 e ciascuna sequenza cosi ottenu- 
ta viene fornita a due terzi filtri FIR di sintesi (e 
20 quindi complessivamente a quattro terzi filtri FIR di 
sintesi H 02 i, H 122 , H 023 , H-- 4 , H 02 5, H 12 €, H 02 i e Hi 23 ) , gene- 
rando otto sequenze di campioni. Le otto sequenze di 
campioni sono quindi fornite a otto terze unita di ri- 
duzione campioni 107-114 ed elaborate in rispettivi 
25 filtri neuro-fuzzy 120-127 aventi la struttura di figu- 
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ra 1. Le sequenze di campioni presenti sulle uscite dei 
filtri neuro-fuzzy 120-127 sono quindi incrementate at- 
traverso unita di incremento campioni 130-137 e fornite 
a rispettivi primi filtri FIR di sintesi G C 2i/ G 12 2, G 02 3 r 
5 G124/ G025/ G126/ G027 e Gi-8- Le sequenze di campioni cosi 
ottenute sono sommate a due a due tramite quattro som- 
matori 140-143 (con un processo inverso a quello subito 
a valle delle seconde unita di riduzione dei campioni 
100-103) , fornite a quattro unita di incremento campio- 

10 ni 146-149 e filtrate nuovamente attraverso quattro . se- 
condi filtri FIR di smtesx Gqu, Gn2; G013 e G114. 

Le sequenze di campioni cosi ottenute sono sommate 
a due a due tramite due sommatori 150, 151 (con un pro- 
cesso inverso a quello subito a valle delle prime unita 

15 di riduzione dei campioni 96, 97), incrementate tramite 
due unita di incremento campioni 154, 155; filtrate at- 
traverso due terzi filtri FIR di sintesi G 0 oi e G102 e in- 
fine sommate in un sommatore 160 in modo da fornire il 
segnale di uscita s 0 (t). 

20 La figura 8 mostra uno schema di flusso relativo 

alia sequenza di passi eseguiti con un filtro multiri- 
soluzione con un numero arbitrario, prefissato, di sud- 
divisioni in sottospazi. 

Inizialmente, passo 160, vengono caricati i cam- 

25 pioni del segnale di ingresso s x (t) ; quindi un contatore 
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di divisione spazi j viene inizializzato a -1, passo 
161, e incremental di una unita, passo 162; un conta- 
tore di sottospazi k viene inizializzato a zero, passo 
163, e incremental di una unita, passo 164. Quindi, i 
5 campioni del segnale di ingresso s x (t) vengono filtrati 
utilizzando il filtro Ho 3k (quindi, alia prima iterazio- 
ne, utilizzando il filtro H„oi) , passo 165, e i campioni 
filtrati vengono ridotti di numero ( "downsampled" ) , 
passo 166. Quindi il contatore di sottospazi k viene 
10 incrementato, passo 170; i campioni del segnale di in- 
gresso Si(t) vengono filtrati utilizzando il filtro H ljk 
(quindi, alia prima iterazione, utilizzando il filtro 
H 1CZ ), passo 171, e i campioni filtrati vengono ridotti 
di numero, passo 172. 
15 Viene verificato se il contatore di sottospazi k e 

pari a 2 l:+1> , in caso negativo viene ripetuto il ciclo 
comprendente i passi 164-172 (nelle iterazioni succes- 
sive alia prima, il filtraggio viene eseguito sui cam- 
pioni ottenuti nella iterazione precedente con il fil- 
20 tro H Cu - 1)(k - 2) o H!,-,,^,), in caso positivo (uscita SI 
dal passo 173) si verifica se la suddivisione in sotto- 
spazi e terminata (cioe se il contatore di divisione 
spazi j ha raggiunto il massimo previsto) , passo 174. 
in caso negativo, si ritorna al passo 162; in caso po- 
25 sitivo, tutte le sequenze di campioni ottenute vengono 
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filtrate utilizzando il filtro neuro-fuzzy di figura 1, 
passo 175 . 

Quindi il contatore di divisione spazi j viene 
inizializzato al suo valore massimo piu una unita, pas- 
5 so 180, e quindi decrementato di una unita, passo 181; 
il contatore di sottospazi k viene inizializzato a ze- 
ro, passo 182, e incrementato di una unita, passo 183. 
Quindi, la prima sequenza di campioni in uscita dal 
primo filtro neuro-fuzzy (120 in figura 7) viene incre- 

10 mentata ( "upsampled" ) , passo 184, e filtrata utilizzan- 
do un filtro Go 3 < (quindi, alia prima iterazione, usandb 
il filtro G 021 ) / passo 185. Quindi il contatore di sotto- 
spazi k viene nuovamente incrementato, passo 190, una, 
seconda sequenza di campioni in uscita da un secondd 

15 filtro neuro-fuzzy (121 in figura 7) viene incremental 
ta, passo 191, e filtrata utilizzando un filtro Gi 3k 
(quindi, alia prima iterazione, usando il filtro G 12 i), 
passo 192. I campioni in uscita dai filtri G 0: k e G-.-, k 
vengono poi sommati, passo 193. 

20 Si verifica quindi se il contatore di sottospazi k 

e pari a 2 (3t11 , passo 194; in caso negativo viene ripe- 
tuto il ciclo comprendente i passi 183-193 (elaborando 
le sequenze di campioni in uscita dai successivi filtri 
neuro-fuzzy) , in caso positivo (uscita SI dal passo 

25 194), si verifica se si e arrivati alia fine, passo 
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195; in caso negativo si ritorna al passo 181, decre- 
mentando il contatore di divisione spazi j ed elaboran- 
do le sequenze di campioni precedentemente incrementa- 
te, filtrate e sommate. L'anello definito dai passi 
5 181-194 viene ripetuto fino ad ottenere una unica se- 
quenza di campioni, corrispondente al segnale di uscita 
s 0 (t), uscita SI dal passo 195. 

La figura 9 mostra 1 ' implementazione hardware del 
filtro neuro-fuzzy 1 di figura 1. In dettaglio, il fil- 
10 tro neuro-fuzzy 1 comprende una memoria dati 200, tre 
memorie di lavoro 201-203, un modulo per il calcolo 
delle caratteristiche di segnale 205, una unita di co- 
mando 206, una unita di memoria delle uscite di primo 
strato 207, una unita di calcolo delle uscite di secon- 
15 do. strato 208, una unita di calcolo dei pesi di rico- 
struzione 209 e un'unita di ricostruzione 210. 

La memoria dati 200 memorizza i (2N+1) campioni 
e(i) di ciascuna finestra di lavoro ed e costituita da 
(2N+1) registri, ciascuno a 16 bit. 
20 Le memorie di lavoro 201-203 sono memorie non vo- 

latili, ad esempio a ROM, PROM, EPROM, EE PROM o flash. 

In particolare, la prima memoria di lavoro 201 me- 
morizza la prima caratteristica di segnale XI (i) e com- 
prende (2N+1) locazioni di memoria a 16 bit. Dato che 
25 il valore della prima caratteristica di segnale XI (i) 
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per il campione i-esimo e costante nelle diverse fine- 
stre di campioni, come risulta evidente dalla (1) , il 
contenuto della prima memoria di lavoro 2 01 non deve 
essere modif icato durante la fase di apprendimento o 
5 durante il f unzionamento del filtro neuro-fuzzy 1 e la 
prima memoria di lavoro 201 puo essere realizzata in 
una qualsiasi delle tecnologie sopra indicate. 

La seconda memoria di lavoro 202 memorizza il va- 
lore delle due curve gaussiane descritte dalla (2) in 

10 funzione del valore delle caratteristiche di segnale 
XI, X2 e X3 . Dato che il valore di tale curve gaussiane 
dipende dai pesi di secondo strato W m (l,k), W v (l r k), nel 
caso sia prevista una fase di apprendimento, la seconda 
memoria di lavoro 202 deve essere di tipo programmabi- 

15 le, • ad esempio di tipo-EPROM, EE PROM o. flash... Per. evi>- 
tare l'uso di una memoria di dimensioni eccessive, le 
funzioni gaussiane (che rappresentano le funzioni di 
appartenenza delle caratteristiche di segnale XI, X2, 
X3, come sopra discusso) sono memorizzate come valori 

20 discreti, in base al livello di accuratezza desiderate 
Ad esempio, nel caso che le funzioni di appartenenza 
assumano 256 valori, con una accuratezza di 16 bit per 
ciascun valore, considerando due gruppi di regole fuzzy 
per ogni caratteristica di segnale XI, X2, X3, la se- 

25 conda memoria di lavoro 202 deve presentare una capaci- 



ta di memoria pari a 256x16x6 bit. La seconda memoria 
di lavoro 202 viene quindi indirizzata a partire dal 
valore corrente (relativo al campione i-esimo) delle 
caratteristica di segnale XI, X2, X3, fornite dalla 
5 prima memoria di lavoro 201 e dal modulo per il calcolo 
delle caratteristiche 205, e fornisce in uscita (al- 
l'unita di memoria delle uscite di primo strato 207) il 
valore delle sei uscite di primo strato oLl (m) . 

La terza memoria di lavoro 203 memorizza i pesi di 
10 secondo strato W FA (m,n) e i pesi di terzo strato W DF (n) . 
Dato che i pesi di terzo strato W DF (n) vengono general- 
mente modificati durante la fase di apprendimento, la 
terza memoria di lavoro 203 e di tipo programmabile, 
analogamente alia seconda memoria di lavoro 202. In 
15 dettaglio, . se si hanno MxN pesi di secondo strato 
W FA (m,n) e N pesi di terzo strato W DF (n), la memoria di 
lavoro 203 comprende (16xMxN+ 16xN) bit. 

II modulo per il calcolo delle caratteristiche 205 
comprende una rete hardware non mostrata, rappresentata 
20 divisa in una prima parte 205a di calcolo della seconda 
caratteristica di segnale X2(i) e in una seconda parte 
205b di calcolo della terza caratteristica di segnale 
X2(i). Nel modulo per il calcolo delle caratteristiche 
205 vengono eseguite le operazioni descritte dalle 
25 equazioni (2) e (3) e comprende una memoria di campioni 



; 29-1 1-2000 

\ 



Rnintedr^e-J 1^20oT| 

_ iasci -art- ;-»T v«l 



ri ■ — i — 

LSPEG, 



00830782J 



■27- 



(avente 2N+1 iocazioni), un'unita per il calcolo del 
valore medio av, un'unita per il calcolo del naassimo 
tra le differenze fra tutti i campioni di ingresso e il 
campione centrale max (diff) , un'unita per il calcolo 
5 del massimo tra le differenze fra tutti i campioni di 
ingresso e il valore medio max (dlff_av) e un'unita per 
il calcolo delle frazioni definite dalle equazioni (2) 
e (3) . 

L'unita di memoria delle uscite di primo strato 

10 207 e formata da sei registri 212 che memorizzano le 
uscite di primo strato oLl (m) fornite dalla seconda me- 
moria di lavoro 202. 

L'unita di calcolo delle uscite di secondo strato 
208 e formata da due moduli, e precisamente un primo 

15 modulo moltiplicativo 208a ed un modulo di ..minimo 208b-. 
In dettaglio, il primo modulo moltiplicativo 208a e 
formato da sei prime unita moltiplicative 213 che mol- 
tiplicano, ciascuna, una rispettiva uscita di primo 
strato oLl (m) (fornita dall'unita di memoria delle 

20 uscite di primo strato 207) con n rispettivi pesi di 
secondo strato W FA (m,n) (forniti dalla seconda memoria 
di lavoro 203); il secondo modulo moltiplicativo 208b e 
formato da sei unita di minimo 214 che, a partire dai 
sei rispettivi prodotti oLl (m) x.W FA (m, n) , ne calcolano il 

25 minimo, fornendo in uscita una rispettiva uscita di se- 



condo strato oL2 (n) . 

L'unita di calcolo dei pesi di ricostruzione 209 e 
formata da due moduli, e precisamente un secondo modulo 
moltiplicativo 209a ed un modulo di defuzzif icazione 
209b. In dettaglio, il secondo modulo moltiplicativo 
209a e formato da sei seconde unita moltiplicative 215 
che moltiplicano una rispettiva uscita di secondo stra- 
to oL2 (n) , fornita dall' unita di calcolo delle uscite 
di secondo strato 208, per un rispettivo peso di terzo 
strato W DF (n), fornito dalla terza memoria di lavoro 
203. Il modulo di defuzzif icazione 209b calcola i pesi 
di ricostruzione oL3, sommando i prodotti forniti dal 
secondo modulo moltiplicativo 209a, sommando fra loro 
le uscite di secondo strato oL2 (n) e f acendo il rappor- 
to fra le due somme, secondo 1'equazione (6). 

L'unita di ricostruzione 210 memorizza i pesi di 
ricostruzione oL3 (i) man mano che questi vengono forni- 
ti dall'unita di calcolo dei pesi di ricostruzione 209 
e, non appena ha immagazzinato 2N+1 pesi di ricostru- 
zione oL3, calcola il campione di uscita u(2N+l), in 
base alia (7), utilizzando anche il valore dei 2N+1 
campioni, forniti dalla memoria dati 200. In seguito, 
al ricevimento del successivo peso di ricostruzione 
oL3(2N+2), esso calcola il campione di uscita u(2N+2), 
usando anche i 2N pesi di ricostruzione ol3 precedenti 
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e altrettanti campioni provenienti dalla memoria dati 
200, in accordo con l'equazione (7) . 

L'unita di controllo 206 determina la sequenza di 
elaborazione e il caricamento/ trasf erimento dei dati 
5 fra le diverse unita e moduli. Per calcolare un singolo 
campione di uscita, l'unita di controllo fa ripetere la 
sequenza dei passi dell ' algoritmo fuzzy per 2N+1 volte, 
aggiorna la memoria dati 200 al termine di 2N+1 cicli e 
controlla il caricamento di successivi 2N+1 campioni. 

10 I vantaggi ottenibili con il metodo e il filtro 

descritto sono i seguenti. In primo luogo, essi consen- 
tono di ridurre il rumore del segnale di ingresso, s'ia 
di tipo bianco, sia di tipo colorato e consentono la 
separazione di segnali con caratteristiche different!. 
.15 II . filtro preserva i front! ripidi- dei. ..segnali ,. senza 
determinare perdite delle caratteristiche del segnale, 
come evidente dal confronto fra le figure 10a, 10b, 
10c. In particolare, la figura 10a e relativa ad un se- 
gnale vocale non rumoroso (e precisamente ad un fram- 

20 mento di segnale corrispondente alia vocale "e", con un 
campionamento a 44,1 KHz e una risoluzione di conver- 
sione A/D di 16 bit) , la figura 10b e relativa alio 
stesso segnale di figura 10a, in presenza di rumore 
bianco, e la figura 10c mostra il risultato del fil- 

25 traggio del segnale di figura 10b con il dispositivo di 



filtro 95 di figura 7. 

Dato che il filtro puo essere addestrato, esso puo 
essere sia adattato ad uno specifico tipo di segnale 
iniziale sia modificato successivamente in caso di ne- 
cessity. Ad esempio, il filtro puo essere adattato ini- 
zialmente ad un primo tipo di segnale acustico (ad 
esempio voce maschile con rumore bianco, eseguendo un 
addestramento con due segnali, uno in presenza di rumo- 
re bianco ed uno in assenza di rumore) . Quindi, la fase 
di apprendimento puo essere ripetuta per segnali di ti- 
po diverso, ad esempio una voce femminile. I nuovi pesi 
cosi calcolati possono sostituire i precedenti, o esse- 
re memorizzati in aggiunta; in questo caso l'unita di 
controllo 206 del filtro 1 puo comandare l'invio alle 
unita 207-209 (figura 9) dei primi o dei secondi cam- 
pioni a seconda dell'uso del filtro. 

Risulta infine chiaro che al metodo e al filtro 
qui descritto ed illustrato possono essere apportate 
numerose modifiche e varianti, tutte rientranti 
nell'ambito del concetto inventivo, come definito nelle 
rivendicazioni allegate. 
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RIVENDICAZIONI 
1. Dispositivo di filtraggio (80; 95) di segnali 
aventi fronti ripidi, caratterizzato dal fatto di com- 
prendere almeno un filtro neuro-fuzzy (1) . 
5 2. Dispositivo di filtraggio secondo la rivendica- 

zione 1, caratterizzato dal fatto che detto filtro neu- 
ro-fuzzy (1) comprende : 

un'unita di calcolo caratteristiche di segnale 
(2), ricevente campioni di ingresso e(i) di un segnale 
10 da filtrare In e generante caratteristiche di segnale 
(XI (i) , X2 (i) , X3 (i) ) ; 

una rete neuro-fuzzy (3) , ■ ricevente dette caratte- 
ristiche di segnale e generante pesi di ricostruzione 
{oL3(i)); e 

15 un'unita di ricostruzione a media mobile' (4), ri- 

cevente detti campioni di ingresso e(i) e detti pesi di 
ricostruzione (oL3(l)) e generante campioni di uscita 
(u(i)) a partire da detti campioni di ingresso e detti 
pesi di ricostruzione. 

20 3. Dispositivo di filtraggio secondo la rivendica- 

zione 2, caratterizzato dal fatto che detta unita di 
calcolo caratteristiche di segnale (2) comprende una 
memoria (200) memorizzante una prima pluralita di cam- 
pioni di ingresso (e(i)) formante una prima finestra ed 

25 una rete di calcolo caratteristiche (201, 205) riceven- 



te detta prima pluralita di campioni di ingresso e for- 
nente una pluralita di caratteristiche (XI (i), X2(i), 
X3(i)) per ogni detto campione di ingresso. 

4. Dispositivo secondo la rivendicazione 3, carat- 
5 terizzato dal fatto che detta rete di calcolo caratte- 
ristiche (201, 205) comprende primi mezzi di determina- 
zione caratteristiche (201) , generanti, per ogni detto 
campione di ingresso (e(i)), una prima caratteristica 
di segnale (XI (i) ) correlata ad una posizione di detto 
.0 campione di ingresso in detta prima finestra; secondi 
mezzi di determinazione caratteristiche (205a), gene- 
ranti, per ogni detto campione di ingresso, una seconda 
caratteristica di segnale (X2(i)) correlata alia diffe- 
renza fra detto campione di ingresso ed un campione 
15 centrale (e (N) ) in detta prima finestra; e terzi mezzi 
di determinazione caratteristiche (205b) , generanti, 
per ogni detto campione di ingresso, una terza caratte- 
ristica di segnale (X3(i)) correlata alia differenza 
fra detto campione di ingresso ed un valore medio di 
20 campione in detta prima finestra. 

5. Dispositivo secondo la rivendicazione 4, carat- 
terizzato dal fatto che detti primi mezzi di determina- 
zione caratteristiche (201) generano detta prima carat- 
teristica di segnale (XI (i)) per un campione di ingres- 
25 so e(i) in base alia relazione: 
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in cui i e la posizione di detto campione di in- 
gresso e(i) in detta prima finestra ed N e la posizione 
di un campione centrale in detta prima finestra. 
5 6. Dispositivo secondo la rivendicazione 5, carat- 

terizzato dal fatto che detti primi mezzi di determina- 
zione caratteristiche comprendono una prima memoria di 
lavoro (201) memorizzante detta prima caratteristica di 
segnale (XI (i) ) per ogni detto campione di ingresso 
10 (e(i)). 

7. Dispositivo secondo una qualsiasi delle riven- 
dicazioni 4-6, caratterizzato dal fatto che detti s.e- 
condi mezzi di determinazione caratteristiche (205a) 
generano detta seconda caratteristica di segnale 

15 (X2(i)) per un campione di ingresso e(i) in base alia 
relazione : 

max(diff) 

in cui e (N) e un campione centrale in detta prima 
finestra e max (diff) e il massimo delle differenze fra 
20 tutti i campioni di ingresso in detta prima finestra e 
detto campione centrale. 

8. Dispositivo una qualsiasi delle rivendicazioni 
4-7, caratterizzato dal fatto che detti terzi mezzi di 
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determinazione caratteristiche (205b) generano detta 

terza caratteristica di segnale (X3(i)) per un campione 

di ingresso e(i) in base alia relazione: 
\e(i)-av\ 
max(diff_av) 

5 in cui av e il valore medio dei campioni di in- 

gresso in detta prima finestra e max (diff_ a v) e il mas- 
simo delle differenze fra tutti i campioni di ingresso 
in detta prima finestra e detto valore medio av. 

9. Dispositivo secondo una qualsiasi delle riven- 
10 dicazioni 2-8, caratterizzato dal fatto che detta rete 
neuro-fuzzy (3) comprende: 

neuroni di f uzzif icazione (6, 202), riceventi der- 
te caratteristiche di segnale (XI (i), X2 (i) , X3(i)) di 
un campione di ingresso (e(i)) e generanti uscite di 
15 primo strato (oil) costituenti il livello di confidenza 
di dette caratteristiche di segnale rispetto a prefis- 
sate funzioni di appartenenza; 

neuroni di fuzzy AND (7) , riceventi dette uscite 
di primo strato e generanti uscite di secondo strato 
20 (oL2) derivanti da regole fuzzy; e 

un neurone di def uzzif icazione (8), ricevente det- 
te uscite di secondo strato e generanti un peso di ri- 
costruzione (oL3) per ogni detto campione di ingresso 
(e(i)) utilizzando il criterio del centro di gravita. 
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10. Dispositivo secondo la rivendicazione 9, ca- 
ratterizzato dal fatto che dette funzioni di apparte- 
nenza sono funzioni gaussiane e dette uscite di primo 
strato sono calcolate in base all ' equazione : 

oL\(l, k) = exp 



r Xl-W m {l,k)^ S 



W v (I,k) 

in cui oLl(l,k) e un'uscita di primo strato; XI e 
una caratteristica di segnale, W m (l,k) e il valore medio 
e W„(l,k) e la varianza di una funzione gaussiana. 

11. Dispositivo di filtraggio secondo la rivendi- 
10 cazione 9 o 10, caratterizzato dal fatto di comprendere 

due funzioni di appartenenza per ogni detta caratteri- 
stica di segnale (XI (i) , X2(i), X3(i)). 

12. Dispositivo di filtraggio secondo una qualsia- 
si delle rivendicazioni 9-11, caratterizzato dal fatto 

15 che detti neuroni di primo strato (6) comprendono una 
seconda raemoria di lavoro (202) memorizzante valori di 
dette uscite di primo strato (oLl(l,k)) per ogni valore 
di dette caratteristiche di segnale (XI (i), X2(i), 
X3(i)). 

20 13. Dispositivo di filtraggio secondo una qualsia- 

si delle rivendicazioni 9-12, caratterizzato dal fatto 
che dette regole fuzzy sono di tipo AND e dette uscite 
di secondo strato (oL2) utilizzano la norma del minimo . 

14. Dispositivo di filtraggio secondo la rivendi- 



- 6 - 



cazione 13, caratterizzato dal fatto che dette uscite 
di secondo strato sono calcolate in base all ' equazione : 
oL2(n) = min On, n) - oL\(m) } 

rt 

in cui oL2(n) e un'uscita di secondo strato, 
W FA (m,n) e un peso di secondo strato e oLl(l,k) e 
un'uscita di primo strato. 

15. Dispositivo di filtraggio secondo la rivendi- 
cazione 14, caratterizzato dal fatto che detti neuroni 
di secondo strato (7) comprendono una pluralita di uni- 
ts moltiplicative (213) ed una pluralita di moduli di 
minimo (214) collegati in cascata. 

16. Dispositivo di filtraggio secondo una qualsia- 
si delle rivendicazioni 9-15, caratterizzato dal fatto 
che detti pesi di ricostruzione sono calcolati in base 

15 all 1 equazione : 

fw DF (n)-oL2(n) 



10 



20 



±oL2(n) 



in cui oL3 e un peso di ricostruzione, W DF (n) sono 
pesi di terzo strato e oL2 (n) e un'uscita di secondo 
strato . 

17. Dispositivo di filtraggio secondo una qualsia- 
si delle rivendicazioni 9-16, caratterizzato dal fatto 
che detta unita di ricostruzione a media mobile (4) ri- 
ceve una seconda pluralita di campioni di ingresso 
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(e(i)) formante una seconda finestra ed una corrispon- 
dente pluralita di pesi di ricostruzione (oL3(i)) e 
calcola ciascun detto campione di uscita (u(i)) in base 
alia equazione: 

f>I3(/-y)-e(/-y) 

7=0 

in cui e(i-j) e un campione di ingresso (i-j)- 
esimo e oL3(±-j) e un peso di ricostruzione associato 
ad un campione di ingresso (i-j)-esimo. 

18. Dispositivo di filtraggio secondo una qualsia- 
si delle rivendicazioni 9-17, caratterizzato dal fatto 
di comprendere una unita di addestramento (5) avente un 
primo ingresso collegato a detta unita di ricostruzione 
a media mobile (4) e ricevente detti campioni di uscita" 
(u(i))/ un secondo ingresso ricevente un segnale di 
uscita desiderato (Tg) ed un'uscita collegata a detta 
rete neuro-fuzzy (3) per fornire valori di pesatura ot- 
timizzati . 

19. Dispositivo di filtraggio secondo una qualsia- 
si delle rivendicazioni 1-18, caratterizzato dal fatto 
di comprendere almeno un primo stadio di divisione 
("splitting stage") (H :( H 0 , 84, 85; H 00 i/ H 102 , 96, 97) 
generante almeno due flussi di campioni da filtrare; un 
detto filtro neuro-fuzzy (86, 87; 120-127) per ogni 



Printed:26-1 1-2001 



iSPEC 



; 00830782 



- 8 - 



flusso di campioni da filtrare, detti filtri neuro- 
fuzzy generando ciascuno un rispettivo flusso di cam- 
pioni filtrati; ed un primo stadio di ricomposizione 
(88, 89, 92, G 0 , Gi; 154, 155, 160, G 0 ci, G 102 ) ricevente 
5 detti flussi di campioni filtrati e generante un unico 
flusso di campioni di uscita (s 0 ) - 

20. Dispositivo di filtraggio secondo la rivendi- 
cazione 19, caratterizzato dal fatto di comprendere 
inoltre una pluralita di ulteriori stadi di divisione 
10 (Hou-Hi28/ 100-114), disposti in cascata fra loro e a 
detto primo stadio di divisione (H 00 i/ Hi 02 , 96, 97) ed 
una pluralita di ulteriori stadi di ricomposizione 
(130-151, G 0 2i-Gii«)/ disposti in cascata fra loro, fra 
detti filtri neuro-fuzzy (120-127) e detto primo stadio 
15 di ricomposizione (154, 155, 160, Gooir Gio?) ■ 

21. Dispositivo di filtraggio secondo la rivendi- 
cazione 19 o 20, caratterizzato dal fatto che detti 
stadi di divisione (H w H Cf 84, 85; Hooi-Hi«, 96-114) com- 
prendono ciascuno un primo ed un secondo filtro di ana- 
20 lisi (H w H 0; Hooi-H 12 3) fra loro in quadratura, riceventi 
un flusso di campioni da dividere e generanti ciascuno 
un rispettivo flusso di campioni divisi e una prima ed 
una seconda unita a riduttore di campioni ( "downsampler 
unit") (84, 85; 96-114) ricevente ciascuna un rispetti- 
25 vo detto flusso di campioni divisi, 
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e dal fatto che detti stadi di ri compos izione (88, 
89, 92, G 0/ G-.; 130-160, G021-G102) comprendono ciascuno 
una prima ed una seconda unita a incremento di campioni 
("upsampler unit") (88, 89; 130-137, 146-149, 154, 
5 155) , riceventi ciascuna un rispettivo flusso di cam- 
pioni da incrementare e generando un rispettivo flusso 
di campioni incrementati; un primo ed un secondo filtro 
di sintesi (G 0 , G a ; G021-G102) fra loro in quadratura e 
complementari a detti filtri di analisi, riceventi cia- 

10 scuno un rispettivo flusso di campioni incrementati ' e 
generando ciascuno un rispettivo flusso di campioni 
parziali; ed un nodo sommatore (92; 140-143, 150, 151, 
160), ricevente detti flussi di campioni parziali e ge- 
nerante un flusso di campioni sommati . 

15 22. Dispositivo di filtraggio secondo la rivendi- 

cazione 21, caratterizzato dal fatto che detti filtri 
di analisi (H*, Ho ; Hooi-Hi^a) sono filtri QMF e detti 
filtri di sintesi (G 0 , Gi; G021-G1C2) sono filtri QMF com- 
plementari a detti filtri di analisi. 

20 23. Dispositivo di filtraggio secondo la rivendi- 

cazione 22, caratterizzato dal fatto che detti filtri 
QMF (Hi, H 0; H001-H128; Go/ Gi; G 0 2i - Gi 0 2) sono filtri convo- 
lutivi . 



24 . 



Metodo per ridurre 



il rumore in 



un filtro do- 



25 



tato di 



fronti ripidi, cara 



tterizzato dal 



fatto di ese- 
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guire una fase di filtrazione neuro-fuzzy. 

25. Metodo secondo la rivendicazione 24, caratte- 
rizzato dal fatto di comprendere le fasi di : 

calcolare caratteristiche di segnale (XI (i), 
5 X2(i), X3(i)) a partire da campioni di ingresso (e(i)) 
di un segnale da filtrare (In); 

calcolare pesi di ricostruzione (oL3) a partire da 
dette caratteristiche di segnale tramite una rete neu- 
ro-fuzzy (3) ; e 

10 ricostruire a media mobile un segnale di uscita 

(Out) formato da una pluralita di campioni di uscita 
(u(i)) a partire da detti campioni di ingresso e detti 
pesi di ricostruzione. 

26. Metodo secondo la rivendicazione 25, caratte- 
15 rizzato dal fatto di comprendere le fasi di : 

memorizzare una prima pluralita di campioni di in- 
gresso (e(i)) formante una prima finestra; 

calcolare una pluralita di caratteristiche di se- 
gnale (XI (i)/ X2(i), X3(i)) per ogni detto campione di 
20 ingresso a partire da detta prima pluralita di campioni 
di ingresso. 

27. Metodo secondo la rivendicazione 26, caratte- 
rizzato dal fatto che detta fase di calcolare una plu- 
ralita di caratteristiche di segnale per ogni detto 

25 campione di ingresso comprende le fasi di : 
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calcolare una prima caratteristica di segnale 
(XI (i) ) correlata ad una posizione di detto campione di 
ingresso (e(i)) in detta prima finestra; 

calcolare una seconda caratteristica di segnale 
(X2(i)) correlata alia differenza fra detto campione di 
ingresso ed un campione centrale in detta prima fine- 
stra; e 

calcolare una terza caratteristica di segnale 
(X3 (i) ) correlata alia differenza fra detto campione di 
ingresso ed un valore medio di campione (av) in detta 
prima finestra. 

28. Metodo secondo la rivendicazione 26 o 27, ca- 
ratterizzato dal fatto che detta fase di calcolare pesi 
di ricostruzione (oL3) comprende le fasi di : 

eseguire un ' operazione di f uzzif icazione e ca'lco*- 
lare uscite di primo strato (oLl) costituenti il livel- 
lo di confidenza di dette caratteristiche di segnale 
(XI (i), X2(i), X3(i)) rispetto a prefissate funzioni di 
appartenenza; 

eseguire un ' operazione di fuzzy AND e generare 
uscite di secondo strato (oL2) derivanti da regole fuz- 
zy a partire da dette uscite di primo strato; e 

eseguire un ' operazione di defuzzif icazione su det- 
te uscite di secondo strato e generare un peso di rico- 
struzione {oL3(l)) per ogni detto campione di ingresso 



(e(i)) utilizzando il criterio del centro di gravita. 

29. Metodo secondo la rivendicazione 28, caratte- 
rizzato dal fatto che dette funzioni di appartenenza 
sono funzioni gaussiane e dette uscite di primo strato 
5 sono calcolate in base all ' equazione : 

'Xl-W m {l,k) 



10 



oL\(J,k) = exp 



in cui oLl(l,k) e un'uscita di primo strato; XI e 
una caratteristica di segnale, W^l.k) e il valore medio 
e W v (1 ,k) e la varianza di una funzione gaussiana. 

30. Metodo secondo la rivendicazione 28 o 29, ca- 
ratterizzato dal fatto che dette uscite di secondo 
strato sono calcolate in base all 1 equazione : 
oL2(«) = min §V FA (m, n) ■ oL\(m)} 

in cui oL2(n) e un'uscita di secondo strato, 
15 W FA (m,n) e un peso di secondo strato e oLl(l,k) e 
un'uscita di primo strato. 

31. Metodo secondo una qualsiasi delle rivendica- 
zioni 28-30, caratterizzato dal fatto che detti pesi di 
ricostruzione sono calcolati in base all ' equazione : 
jy DF (n)-oL2(ri) 



20 oL3=— 



in cui oL3 e un peso di ricostruzione, W DF (n) sono 
pesi di terzo strato e oL2 (n) e un'uscita di secondo 
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32. Metodo secondo una qualsiasi delle rivendica- 
zioni 26-31, caratterizzato dal fatto che detta fase di 
ricostruire a media mobile comprende le fasi di : 

ricevere una seconda pluralita di campioni di in- 
gresso (e(i)) formante una seconda finestra ed una cor- 
rispondente pluralita di pesi di ricostruzione {oL3(±)) 
e calcolare ciascun detto campione di uscita in base 
alia equazione: 

f>Z3(/-j) •*('-./') 

s-\ J=° 

"(0 = 5iv 

Z *(<-./) 

in cui-e(i-j) e campione di ingresso ( i- j ) -esimd; e 
oL3(i-j) e un peso di ricostruzione associato ad un 
campione di ingresso (i-j)-esimo. 

33. Metodo secondo una qualsiasi delle rivendica- 
zioni 24-32, caratterizzato dal fatto di effettuare una 
fase di addestramento di pesi utilizzati in detta fase 
di filtrazione neuro-fuzzy. 

34. Metodo secondo la rivendicazione 33, caratte- 
rizzato dal fatto che detta fase di addestramento com- 
prende le fasi di : 

generare un segnale di ingresso di conf igurazione 

no t a ; 

filtrare detto segnale di ingresso di configura- 
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zione nota per ottenere un segnale di uscita di prova; 

confrontare detto segnale di uscita di prova con 
un segnale desiderato (Tg) per ottenere una distanza 
fra detto segnale di uscita di prova e un segnale desi- 
5 derato; 

calcolare una funzione di conformita ("fitness") a 
partire da detta distanza; 

ottimizzare detti pesi in base a detta funzione di 
conformita . 

1 0 35. Metodo secondo una gualsiasi delle rivendica- 

zioni 26-34, caratterizzato dal fatto di effettuare una 
analisi multirisoluzione con suddivisione di segnale in 
sottobande tramite ondine ortonormali. 

36. Metodo secondo la rivendicazione 35, caratte- 
15 rizzato dal fatto di comprendere le fasi di : 

dividere un flusso di campioni di ingresso in al- 
meno due flussi di campioni da filtrare; 

filtrare ogni flusso di campioni da filtrare tra- 
mite un rispettivo filtro neuro-fuzzy (86, 87; 120-127) 
20 per ottenere almeno due flussi di campioni filtrati; e 

ricomporre detti flussi di campioni filtrati per 
generare un unico flusso di campioni di uscita (s 0 ) . 

37. Metodo secondo la rivendicazione 36, caratte- 
rizzato dal fatto che, prima di eseguire detta fase di 
25 filtrare, detta fase di dividere viene ripetuta un nu- 
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mero prefissato di volte e che, dopo aver eseguito det- 
ta fase di filtrare, detta fase di ricomporre viene ri- 
petuta detto numero prefissato di volte. 

38. Metodo secondo la rivendicazione 36 o 37, ca- 
ratterizzato dal fatto che detta fase di dividere com- 
prende le fasi di : 

alimentare un flusso di campioni da dividere a due 
filtri di analisi fra loro in quadratura; 

generare due flussi di campioni divisi filtrati; e 

ridurre in campioni ( "downsampling" ) detti flussi 
di campioni divisi filtrati, 

e dal fatto che detta fase di ricomporre comprende 
le fasi di : 

incrementare flussi di campioni da incrementare, 
generando flussi di campioni incrementati; 

filtrare detti flussi di campioni incrementati 
tramite due filtri di sintesi fra loro in quadratura e 
complementari a detti filtri di analisi, generando 
flussi di campioni parziali; e 

sommare coppie. di flussi di campioni parziali e 
generare un flusso di campioni sommati. 



RIASSUNTO 

II dispositivo di filtro (80) comprende un filtro neu- 
ro-fuzzy (1; 80) ed implementa una tecnica di filtrag- 
gio a media mobile in cui i pesi per la ricostruzione 
5 finale del segnale (oL3(i)) sono calcolati in una rete 
neuro-fuzzy (3) in base a specifiche regole fuzzy. Le 
regole fuzzy operano su tre caratteristiche di segnale 
(XI (i), X2(i), X3(i)) per ogni campione di ingresso 
(e(i)). Le caratteristiche di segnale sono correlate 
10 alia posizione del campione nella finestra di campioni 
considerata, alia differenza fra un campione e il cam- 
pione al centro della finestra e alia differenza tra un 
campione e la media di campioni nella finestra. II di- 
spositivo di filtro per l'analisi di un segnale vocale 
15 comprende un banco di filtri neuro-fuzzy (86, 87) . II 
segnale viene separato in diverse sotto-bande, secondo 
la teoria delle ondine, utilizzando un banco di filtri 
di analisi, includente una coppia di filtri FIR QMF (H 0 , 
H a ) ed una coppia di riduttori di campioni (85, 86); 
20 ciascun segnale di sotto-banda viene filtrato da un 
filtro neuro-fuzzy (86, 87) e quindi le varie sotto- 
bande vengono ricostruite da un banco di filtri di sin- 
tesi, includente una coppia di incrementatori di cam- 
pioni (88, 89), una coppia di filtri FIR QMF (G 0 , Gi) ed 
25 un nodo sommatore (92) . Figure 1 e 4 
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FILTERING DEVICE AND METHOD FOR REDUCING NOISE IN ELECTRICAL 
SIGNALS, IN PARTICULAR ACOUSTIC SIGNALS AND IMAGES 

The present invention refers to a filtering device and 
5 method for reducing noise in electrical signals, in 

particular acoustic (voice) signals and images. 

As is known, filters operating according to various linear 
and nonlinear techniques are used to remove undesired 
10 components from electrical signals. In particular, undesired 

components may be any type of noise (white noise, flicker, 
etc.) or other types of superimposed acoustic or visual 
signals . 

15 Linear filters are at present the most widely used 

instruments for filtering noise. Finite-impulse filters 
(FIRs) eliminate all the harmonics of a signal having a 
frequency higher than the cutoff frequency of the filter and 
improve the signal-to-noise ratio (SNR) . Another linear 

20 filtering technique is based on the fast Fourier transform 

(FFT) , where the signal is transformed into the frequency 
domain, the undesired harmonics are removed, and then the 
inverse Fourier transform is calculated. 

25 As far as nonlinear techniques are concerned, average 

filters are simple to design and may be implemented through 
simple hardware circuits. Average filters are based on the 
comparison of the individual signal samples in an time 
interval with the average of all the samples in the same 

30 time interval. On the basis of this comparison, the 

individual samples are selectively attenuated. 



35 



All these methods share the disadvantage that, when removing 
the noise, also some of the components of the original 
signal are removed. 
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Furthermore, none of the current techniques, whether linear 
or nonlinear ones, including average filtering, is able to 
preserve steep edges of the signal. If a moving-average 
filter is used, the width of the window must be very small 
if steep edges are to be preserved. However, if the size of 
the window becomes small, there is no significant reduction 
in noise energy. If linear filters are used, all the 
frequencies above the cutoff frequency are eliminated, with 
consequent marked distortion of the signal. 

The aim of the invention is thus to provide a filtering 
method and device that does not cause a sensible 
deterioration of the signal and at the same time preserves 
the edges of the signal. 

According to the present invention a device and method are 
provided for filtering electrical signals, as respectively 
defined in Claim 1 and Claim 24. 

The filter and the method described are particularly useful 
in the case of signals having steep edges, for which the aim 
is to preserve the edges of the signals. Furthermore, it is 
possible to filter signals affected by white and non-white 
noise, such as flicker noise. Through the present method it 
is moreover possible to eliminate from a signal other 
signals that are superimposed on it and are characterized by 
a wide spectral range. 

The device and the method described are based upon a neuro- 
fuzzy network. They are implemented with a moving-average 
filtering technique in which the weighting factors (or 
weights) for the final reconstruction of the signal are 
calculated in a neuro-fuzzy network according to specific 
fuzzy rules. This enables a better reduction of the noise. 
The fuzzy rules operate on different variables, referred to 
as signal features. Described hereinafter are three signal 
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features and six fuzzy rules. 



The proposed filter is suitable for visual signals or 
acoustic signals, even ones with sudden variations. Various 
types of functions or signal features can be used to create 
the rules. With the method described, the signal features 
are correlated to the position of the sample in the 
considered sample window, to the difference between a given 
sample and the sample at the center of the window, and to 
the difference between a given sample and the average of 
samples in the window. These signal features may have a 
considerable influence on the values of the weights for the 
reconstruction of the signal; in addition, they may be 
calculated in a relatively simple manner. 

The method and the filter according to the invention 
moreover comprise a neuro-fuzzy filter bank. In this way, 
the signal may be split into different sub-bands according 
to wavelet theory: each sub-band signal may be filtered by. a 
neuro-fuzzy network, and then the various sub-bands can be 
reconstructed by the synthesis filter bank. As is known from 
wavelet theory, in the first sub-band the signal features 
have a low frequency, whereas in the last sub-band the 
signal features have the maximum frequency. If non-white 
noise (colored noise) is to be removed, this is approximated 
by white noise in each individual sub-band. Given that a 
neuro-fuzzy network works well on white noise, this solution 
leads to a sensible reduction in noise. 

The network is trained by supplying some configurations of 
input and output signals {the configuration of the output 
signal that it is required to obtain as a result of the 
network evolution is called target configuration) . The 
training algorithm is based upon one of the known learning 
methods, such as gradient descent, a genetic algorithm, the 
simulated annealing method, random search, or any other 
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method for function optimization. 

For an understanding of the present invention, preferred 
embodiments thereof are now described, purely to furnish 
5 non-limiting examples, with reference to the attached 

drawings, in which: 

Figure 1 is a block diagram showing the general 
architecture of the filter according to the invention; 

- Figure 2 represents the topology of a block of Figure 1 
10 for a neuro-fuzzy network; 

- Figures 3a and 3b show a flowchart of the operations 
performed by the filter of Figure 1; 

- Figure 4 shows a block diagram of a filtering device using 

the filter of Figure 1; 
15 - Figure 5 shows the transfer functions of digital filters 

belonging to the filtering device of Figure 4; 

- Figure 6 shows the block diagram of some digital filters 
belonging to the filtering device of Figure 4; 

- Figure 7 shows a block diagram of another filtering device 
20 using the filter of Figure 1; 

- Figure 8 shows a flowchart of the operation of the 
filtering device of Figure 7; 

- Figure 9 shows a block diagram for the hardware 
implementation of the filter of Figure 1; and 

25 - Figures 10a, 10b, and 10c respectively show a voice signal 

free from noise, the same signal superimposed on white 
noise, and the same signal after filtering according to the 
invention. 

30 Figure 1 shows a filter 1 comprising a signal-feature 

computation unit 2, a neuro-fuzzy network 3, a 
reconstruction unit 4, and a training unit 5. 



35 



The signal-feature computation unit 2 receives at the input 
a signal In including of a plurality of input samples e(i), 
stores, at each clock cycle, (2N + 1) input samples e(i) 
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(which represent a work window for filter 1) in an internal 
buffer, computes the signal features XI (i), X2(i), and X3 (i) 
for each input sample e(i) on the basis of all the input 
samples belonging to the work window (as described in detail 
5 hereinafter), and supplies the signal features XI (i), 

X2(i), and X3 (i) thus calculated to the neuro-fuzzy network 
3. 



The neuro-fuzzy network 3 processes the signal features 
10 XI (i), X2(i), and X3(i) and generates at the output a 

reconstruction weight oL3(i) for each input sample e(i) . To 
this aim, starting from the signal features XI (i) , X2 (i) , 
and X3(i) and for each input sample e(i), the neuro-fuzzy 
network 3 first performs a fuzzif ication operation, then 
15 applies preset fuzzy rules, and finally carries out a 

defuzzif ication operation. The reconstruction weight oL3 (i) 
thus obtained is hence the weighted sum of all the input 
samples e(i) in the same work window, as explained in detail 
hereinafter with reference to Figures 2 and 3. 

2 0 

The reconstruction unit 4 receives the reconstruction 
weights oL3{i) and the input samples e(i) and, after 
accumulating a sufficient number of input samples e(i) and 
of corresponding reconstruction weights oL3(i), generates an 
25 output sample u(i) the sequence of which forms an output 

signal Out. 

The training unit 5 is operative only initially, so as to 
train the neuro-fuzzy network 3 and modify the weights of 
30 the network with the aim of obtaining an optimal behavior of 

the filter 1, as described in detail hereinafter. 

The signal features XI (i) , X2 (i) , and X3 (i) , computed in the 
signal-feature computation unit 2, are correlated, 
35 respectively, to the distance between each sample and the 

central sample of the considered window, to the difference 
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between a given sample and the sample at the center of the 
window, and to the difference between a given sample and the 
average of samples in the window, and are normalized so as 
to obtain values between 0 and 1. 

5 

In detail, given a window of (2N + 1) input samples e(i), 
with i" = 0, . .., 2N, the signal features XI (i) , X2(i), and 
X3{i) for each input sample e(i) are defined as 

\i-N\ m 

10 X2(0 " max(^) 

_ (3) 
max(<2?# r _av) 

where N is the position of a central sample e(N) in the work 
window; 

max (diff) = max (e(Jc)-e(W)) with Jc = 0, 2N, i.e., the 

15 maximum of the differences between all the input samples 

e(Jt) and the central sample e(N); 

av is the average value of the input samples e(i); and 
max (d±ff_av) = max (e(k)-av) withk= 0, 2N, i.e., the 

maximum of the differences between all the input samples 
20 e{k) and the average value av. 

The neuro-fuzzy network 3 is a three-layer fuzzy network the 
functional representation of which appears in Figure 2, in 
which, for reasons of simplicity, the index i in parenthesis 
25 for the specific sample within the respective work window is 

not indicated. Nevertheless, as- mentioned above, and as will 
emerge more clearly from the flowchart of Figures 3a, 3b, 
the neuro-fuzzy processing represented by Figure 2 is 
repeated for each input sample e(i). 



30 



In detail, starting from the three signal features XI, X2 
and X3 (or generically from 2 signal features XI) and given 
k membership functions of a gaussian type for each signal 
feature (described by the average value W m (l, k) and by the 
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variance W v (l, k) , a fuzzif ication operation is performed in 
which the membership level of the signal features XI, X2 and 
X3 is evaluated with respect to each membership function 
(here two for each signal feature, so that k = 2 ; in all 
5 there are M = 1 x k = 6 membership functions) . 

In Figure 2, the above operation is represented by six 
first-layer neurons 6, each of which, starting from the 
three signal features XI, X2 and X3 (generically designated 
10 by XI) and using as weights the average value W m (l, k) and 

the variance W V {1, k) of the membership functions, supplies a 
first-layer output oil (1 r k) (hereinafter also designated by 
oLl (m) ) calculated as follows : 



15 



30 



oLl(l, k) — oL\(m) - exp 



W v (l,k) 



(4) 



Hereinafter, a fuzzy AND operation is performed, using the 
norm of the minimum, in such a way as to obtain N second- 
layer outputs oL2 (n) . For example, N is equal to 6. As is 
known, a fuzzy AND operation using the norm of the minimum 
20 is based upon rules of the type: 



if X™ is 4° and X™ is A™ and X™ is A™ then 0 is 



if X^ is A™ and AT<"> is and Jf 3 (B) is A™ then 0 is 

in which A^ , A^ B[ X) ,, etc. are linguistic terms, such as 

25 "high" and "low", and the value of the output O for each 

rule is given by the minimum of the membership levels 
multiplied by a weight. 



In practice, with the neuro-fuzzy network of Figure 2, each 
second-layer output oL2(n) is equal to the minimum one among 
the products of the M outputs oLl {m) of the first-layer 
neurons 6 and a respective second-layer weight W FA (m, n) . 



35 



In Figure 2 the above operation is represented by N second- 
layer neurons 7 which implement the equation 
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oL2(«) = min{^(m, /?) -oIA(m)} ( 5 ) 

n 

Finally, the third layer corresponds to a def uzzif ication 
operation and supplies at output a discrete-type 
reconstruction weight oL3, using N third-layer weights 
W DF {n) . The defuzzif ication method is that of the center of 
gravity (centroid) and is represented in Figure 2 by a 
third-layer neuron 8 supplying the reconstruction weight oL3 
according to the equation 



YW DF {n)-oL2{n) 
f>L2(n) 



ra «=i (6) 

10 oL3= ^ 



The reconstruction unit 4 then awaits a sufficient number of 
samples e(i) and corresponding reconstruction weights oL3(i) 
(at least 2N + 1, corresponding to the width of a work 
15 window) and calculates an output sample u(i) as the weighted 

sum of the input samples e(i - j) , with j = 0 ... 2N, using 
the reconstruction weights oL3(i-j) f according to the 
equation 



^oL3(i-j)-e(i-j) 

!>(/-./) 



;=0 



The training unit 5 operates only in an initial learning 
stage, when an input signal In having a known configuration 
is supplied to the filter 1, the output signal Out obtained 
is compared with a target signal Tg, and the distance 
25 between the obtained signal Out and the target signal Tg is 

evaluated on the basis of a fitness function. This fitness 
function may be, for example, the quadratic signal/noise 
ratio having the following expression: 

in which T is the total number of input samples e(i). 



8 
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Using the fitness function and applying a method for 
function optimization, such as the gradient-descent method, 
a genetic algorithm, the simulated- annealing method, and 
5 random search, the first-layer weights (mean value W m {l, k) 

and variance W V {1, k) of the gaussian membership functions) 
and the third-layer weights W DF (n) are modified, and a new 
fitness evaluation is performed. The second-layer weights 
W F a{bi, n) are instead randomly initialized and are not 
10 modified. The learning process is iterated until a preset 

value of the fitness function is achieved or until a preset 
number of genetic algorithms have been generated or a preset 
number of steps of the selected optimization algorithm has 
been performed. 

15 

In this way, the neuro-fuzzy network 3 implements an 
adaptive algorithm and overcomes the limitations of neural 
networks or of fuzzy systems considered separately. In fact, 
fuzzy systems do not have learning capabilities and, if the 
20 selection of the fuzzy rules is not accurate, the fuzzy 

algorithm does not behave satisfactorily. Instead, using the 
neuro-fuzzy network 3 and carrying out a prior learning step 
it is possible to approach a signal having a complex 
mathematical structure, for example the voice in a noisy 
25 environment, without any prior knowledge of the mathematical 

laws governing the said system. 

Operation of the filter 1 of Figure 1 is described in detail 
hereinafter with reference to Figures 3a and 3b. 

30 

Initially, the filter 1 is initialized, step 10, as is a 
window counter p, step 11. Then the window counter p is 
incremented by (2N + 1), equal to the number of input 
samples e(i) in a work window, step 12, and the input 
35 samples of a work window ranging between e(p) and e (p + 2N) 

are loaded, step 13. A sum-of-samples counter h is 



1 1 
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initialized at 0, and a sum-of-samples variable S(0) is 
initialized with the value of the input sample e(p>, step 
14. Next, the sum-of-samples counter h is incremented by one 
unit, step 15, and the sum-of-samples variable S(h) is 
incremented with the value of a sample e (p + h) , step 20. 
After the sum of (2N + 1) input samples (output YES from 
step 21), the sum-of-samples variable S (h) is equal to the 
sum of all the input samples e(i) and may be directly used 
for calculating the signal features. 



Subsequently, using equations (1), (2) and (3) and the value 
of the sum-of-samples variable S (h) , the signal features XI, 
X2, X3 are calculated for each sample e (p) - e (p + 2N) 
belonging to the considered work window, step 22. A sample 
15 counter i for the input samples within the considered work 

window is then reset to zero, step 23, and subsequently 
incremented by one unit, step 24; a feature counter 1 for 
the number of signal features used (in the present case, 3) 
is reset to zero, step 25, and subsequently incremented, 
step 26; and a function counter k for the number of 
membership functions used for each signal feature (in the 
present case, 2) is reset to zero, step 30, and subsequently 
incremented, step 31. Next, the first-layer weights (the 
mean value W m U, k) and the variance W V {1, *> ) are loaded, 
step 32, and the first-layer outputs oil a, k) are 
calculated in accordance with equation (4), step 33. 

Loading of the first-layer weights W m {l, k) , W v [l, k) and 

calculation of the first-layer outputs oLl (1, k) are 

repeated for all the membership functions (output YES from 

step 34) and for all the signal features (output YES from 
step 35) . 



Next, a second-layer neuron counter n is reset to zero, step 
39, and incremented by one unit, step 40; an algorithm 
counter m is reset to zero, step 41, and incremented, step 
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42; and the second-layer weights W FA (m, n) for the n-th 
second-layer neuron 7 are loaded, steps 43, 44. Then, the 
second-layer output oI,2(n) for the n-th second-layer neuron 
7 is calculated using equation (5) , step 45, and the 
5 procedure of steps 40-45 is repeated for all the n second- 

layer neurons 7 . 

At the end (output YES from step 50) , the second-layer 
neuron counter n is reset again, step 51, and incremented by 
10 one unit, step 52; all the third-layer weights W DF {n) are 

loaded in succession, steps 53, 54; at the end of loading 
(output YES from step 54) the reconstruction weight oL3(i) 
is calculated in accordance with equation (6), step 55. 

15 The reconstruction weights o£3(i) thus calculated and the 

respective input samples e(i) are stored, step 60, and the 
procedure described for steps 24-60 is repeated for 
successive input samples e(i) until (2W + 1) input samples 
e(i) are reached (output YES from step 61) . [ 

20 V 

Next, a sum counter j is reset to zero, step 62, and the 
input sample e(i) is added to the previous 2N input samples 
e(i - j) , and the sum is stored in a variable E, step 63. 

25 After the sum of (2N + 1) input samples has been carried 

out, an output sample u(i) is calculated in accordance with 
equation 7, step 65. 

The entire cycle described by steps 12-65 is repeated as 
30 long as input samples e(i) are present; at the end (output 

YES from step 70) , the sample processing procedure 
terminates . 

According to another aspect of the invention, filtering is 
35 based on a multi-resolution analysis obtained through a 

filter bank in phase quadrature. Wavelet theory furnishes 
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the theoretical basis for multi-resolution. 

As is known, a multi-resolution analysis defines a set of 
nested subspaces of a square summable function space, i.e., 

5 the space of the finite-energy functions, widely known in 

physics and electrical engineering. On the basis of the 
above analysis, a projection of a function in one of these 
subspaces isolates the -roughest" details of the function, 
whilst projection of the function in the orthonormal 

10 complement of this subspace isolates the -finer" details of 

the function. The procedure may be iterated to obtain a 
pyramid. From wavelet theory it is known that the operation 
may be performed using a bank of FIR filters, in which each 
FIR filter is followed by a subsampler) . The signal thus 

15 split into sub-bands can be reconstructed using a bank of 

complementary filters, each of which is provided with a 
sample incrementer (upsampler) . 



A first embodiment of the above-mentioned solution is shown 
in Figure 4 and regards a multi-resolution filter 80 
generating two subspaces and using a first pair of Finite 
Impulse Response Quadrature Mirror Filters (FIR QMFs) for 
the analysis and a second pair of FIR QMFs for the 
synthesis . 



in detail, the multi-resolution filter 80 receives, on an 
input 81, an input signal Si (t). The input signal Si (t) is 
supplied to two input FIR filters fib, Hi which perform a 
convolution operation so to output a corresponding filtered 
30 signal e 0 , e, equal to: 



r = 0 
M 



e,(t)=j:M^-Kt-r) (10) 



r=0 



in which M is the order of the filters H 0l Hi; e 0 (t), e x (t) 
is the t-th sample of the respective output sequence; Si(t) 
35 is the t-th sample of the input sequence; h 0 (r) , h 2 (x) is the 



t-th tap of the input FIR filter Ho, Hi, in which 
h: (r) = (-1) h 0 (2 J - r + 1) 
where J is an integer. 

The input FIR filters H Q , H x have transfer functions shown in 
Figure 5 and thus define, respectively, a low-pass filter 
and a " high-pass filter, and have the structure shown in 
Figure 6, where the generic tap h{r) corresponds to h 0 {r) or 
to hi(r) , according to whether it is the input FIR filter H 0 
or the input FIR filter H x . 

The outputs of the input FIR filters H 0 , Hl are each 
connected to a respective subsampling unit 84, 85 which 
discards the odd samples from, the input signal e 0 (£), ej (t) 
and retains only the even samples, generating a respective 
signal e 0 (t), e x (t). The outputs of the subsampling units 84, 
85 are each connected to a respective neuro-fuzzy filter 86, 

87. Both of the neuro-fuzzy filters 86, 87 have the 
structure shown in Figure 1. 

The output signals u 0 (t), Ui(t) of the neuro-fuzzy filters 
86, 87 are each supplied to a respective upsampler 88, 89 
which generates a respective output signal Uou(t), Ui u (t) by 
entering a zero sample between each pair of samples of the 
respective output signal u 0 (t), Ui(t) of the neuro-fuzzy 
filters 86, 87. The outputs of the of the upsampling units 

88, 89 are each connected to a respective output FIR filter 
Go, Gi. These filters too have each a respective transfer 
function given by equation (9) and equation 
(10) , respectively. 

Finally, the output signals of the output FIR filters Go, G x 
are added together, sample by sample, by an adder 92. 

Figure 7 shows a multi-resolution filter 95 using eight 
subspaces. In detail, the input signal Sj.(t) is initially 
supplied to two first synthesis FIR filters H 0 oi/ H102/ 



respectively of the low-pass type and of the high-pass type, 
and is then subsampled in two first subsampling units 96, 
97, in a similar way as described for the units 85, 86 of 
Figure 4 . The sequences of samples thus obtained are each 
supplied to two synthesis filters (and hence altogether to 
four second synthesis FIR filters Hon, H112/ Hoi 3/ and H14) . 
The outputs of the second synthesis FIR filters Hon, Hn 2 , 
JT 013/ and Hu« are then supplied to four second subsampling 
units 100-103, and each sequence thus obtained is supplied 
to two third synthesis FIR filters (and hence altogether to 
eight third synthesis FIR filters H 02 i, H i22 , H 023 , H 124 , H 02 5, 
H126, H 02 7, Hub), generating eight sequences of samples. The 
eight sample sequences are then supplied to eight third 
subsampling units 107-114 and processed in respective neuro- 
fuzzy filters 120-127 having the structure illustrated in 
Figure 1. The sample sequences present on the outputs of the 
neuro-fuzzy filters 120-127 are then incremented via 
upsampling units 130-137 and supplied to respective first 
synthesis FIR filters G 02 i, G i22 , Go 2 3, Gi 24 , G C2 s/ Gi 26 , G 02 7, and 
G 12 8- The sample sequences thus obtained are added up two by 
two through four adders 140-143 (with a reverse process with 
respect to the one followed downstream of the second 
subsampling units 100-103) , supplied to four upsampling 
units 146-149, and filtered again through four second 
synthesis FIR filters G 0 n, G112, G013/ and Gn 4 . 

The sample sequences thus obtained are added up two by two 
through two adders 150, 151 (according to a reverse process 
with respect to the one followed downstream of the first 
subsampling units 96, 97), incremented by two upsampling 
units 154, 155, filtered through two third synthesis FIR 
filters G001 and G 102 , and finally summed in an adder 160 so 
as to supply the output signal s c (t) . 

Figure 8 shows a flowchart of the sequence of steps 
performed using a multi-resolution filter with a preset 



arbitrary number of splittings into subspaces . 

Initially, the samples of the input signal Si(t) are loaded, 
step 160; then, a space split counter j is initialized at — 1 , 
step 161, and incremented by one unit, step 162; a subspace 
counter k is initialized at zero, step 163, and incremented 
by one unit, step 164. Then the samples of the input signal 
Si(t) are filtered using the filter H 0 jk (thus, at the first 
iteration, using the filter ffooi) i step 165, and the filtered 
samples are downsampled, step 166. Next, the subspace 
counter k is incremented, step 170; the samples of the input 
signal Si(t) are filtered using the filter Hi JJt (thus, at the 
first iteration, using the filter H102) / step 171, and the 
filtered samples are downsampled, step 172. 

It is checked whether the subspace counter k is egual to 
2 (J+1] ; if it is not, the cycle comprising steps 164-172 is 
repeated (in the iterations following on the first, 
filtering is performed on the samples obtained in the 
previous iteration with the filter H 0 (j-i! u-2> or Hi (j _d u _ 2 ) ) ; if 
it is (output YES from step 173), it is checked whether .the 
splitting into subspaces is terminated (i.e., whether the 
space split counter j has reached the preset maximum) , step 
174. If it has not, the procedure returns to step 162; if it 
has, all the obtained sample sequences are filtered using 
the neuro-fuzzy filter of Figure 1, step 175. 

Next, the space split counter j is initialized at its 
maximum value plus one unit, step 180, and then decreased by 
one unit, step 181, the subspace counter k is initialized at 
zero, step 182, and incremented by one unit, step 183. Next, 
the first sample sequence at output of the first neuro-fuzzy 
filter (120, in Figure 7) is upsampled, step 184, and 
filtered using a filter G 0jk (thus, at the first iteration, 
using the filter G 02 \) , step 185. Then, the subspace counter 
k is incremented again, step 190; a second sample sequence 
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at the output of a second neuro-fuzzy filter (121, in Figure 
7) is upsampled, step 191, and filtered using a filter G ljk 
(thus, at the first iteration, using the filter G121) / step 
192. The samples at the output of the filters G oj k and G ljk 
5 are then summed, step 193. 

It is then checked whether the subspace counter k is equal 
to 2 t ' +1) , step 194; if it is not, the cycle comprising steps 
183-193 is repeated (processing the sample sequences at the 

10 output of the subsequent neuro-fuzzy filters) ; if it is 

(output YES from step 194), it is checked whether the end 
has been reached, step 195; if it has not, the procedure 
returns to step 181, decreasing the space split counter j 
and processing the sample sequences previously upsampled, 

15 filtered and summed. The loop defined by steps 181-194 is 

repeated until a single sequence of samples is obtained, 
corresponding to the output signal s 0 (t), output YES from 
step 195. 

20 Figure 9 shows the hardware implementation of the neuro- 

fuzzy filter 1 of Figure 1. In detail, the neuro-fuzzy 
filter 1 comprises a data memory 200, three work memories 
201-203, a signal features calculation module 205, a control 
unit 206, a first-layer output memory unit 207, a second- 

25 layer output calculating unit 208, a reconstruction-weight 

calculating unit 209, and a reconstruction unit 210. 

The data memory 200 stores the (2M + 1) samples e(i) of each 
work window and comprises (2N + 1) registers, each having 16 
30 bits. 

The work memories 201-203 are nonvolatile memories, for 
example ROM, PROM, EPROM, EEPROM or flash memories. 

35 In particular, the first work memory 201 stores the first 

signal feature XI (i) and comprises (2N + 1) sixteen-bit 
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memory locations. Since the value of the first signal 
feature XI (i) for the i-th sample is constant in the various 
sample windows, as is evident from equation (1)/ the 
contents of the first work memory 201 must not be modified 
during the learning step or during operation of the neuro- 
fuzzy filter 1, and the first work memory 201 may be 
implemented using any one of the technologies referred to 
above . 

The second work memory 202 stores the values of the two 
gaussian curves described by equation (2) according to the 
values of the signal features XI, X2, and X3 . Since the 
values of these gaussian curves depend upon the second-layer 
weights W m {l, k) , W v (l, k) , when a learning step is provided, 
the second work memory 202 must be of the programmable type, 
for example of the EPROM, EE PROM or flash type. To avoid the 
use of a memory of excessive size, the gaussian functions 
(which represent the membership functions of the signal 
features XI, X2 , X3 , as discussed above) are stored ,as 
discrete values, according to the desired level of accuracy. 
For example, if the membership functions have 256 values, 
with an accuracy of 16 bits per value, considering two sets 
of fuzzy rules for each signal feature XI, X2, and X3, the 
second work memory 202 must have a storage capacity of 
256 x 16 x 6 bits. The second work memory 202 is then 
addressed starting from the current values (corresponding to 
the i-th sample) of the signal features XI, X2, X3 supplied 
by the first work memory 201 and by the signal features 
calculation module 205, and outputs (to the first-layer 
output memory unit 207) the values of the six first-layer 
outputs oLl (m) . 

The third work memory 203 stores the second-layer weights 
W FA (m, n) and the third-layer weights W DF (n) . Since the third- 
layer weights W DF (n) are generally modified during the 
learning step, the third work memory 203 is of a 
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programmable type, as is the second work memory 202. In 
detail, if there are M x N second-layer weights W fA (m, n) and 
N third-layer weights W„ e {n) , the work memory 203 comprises 
(16 x M x N + 16 x N) bits. 

The signal features calculation module 205 comprises a 
hardware network (not shown) represented as divided into a 
first part 205a for calculating the second signal feature 
X2(i) and a second part 205b for calculating the third 
signal feature X2 (i) . In the signal features calculation 
module 205 the operations represented by equations (2) and 
(3) are performed, and this module comprises a sample memory 
(having 2N + 1 locations) , a unit for calculating the 
average value av, a unit for calculating the maximum of the 
differences between all the input samples and the central 
sample max (diff) , a unit for calculating the maximum of the 
differences between all the input samples and the average 
value max(dix--f__av) , and a unit for calculating the fractions 
defined by equations (2) and (3) . 

The first-layer output memory unit 207 comprises six 
registers 212 which store the first-layer outputs oLl (m) 
supplied by the second work memory 2 02. 

The second-layer output calculating unit 208 comprises two 
modules, namely, a first multiplication module 208a and a 
minimum module 208b. In detail, the first multiplication 
module 208a includes six first multiplication units 213 each 
of which multiplies a respective first-layer output oLl (m) 
(supplied by the first-layer output memory unit 207) by n 
respective second-layer weights W FA (m, n) (supplied by the 
second work memory 203); the second multiplication module 
208b includes six minimum units 214 which, starting from the 
six respective products oLl (m) x W FA (m, n) , calculate the 
minimum thereof, supplying at the output a respective 
second-layer output oL2(n). 



The reconstruction-weight calculating unit 209 comprises two 
modules, namely, a second multiplication module 2 09a and a 
defuzzif ication module 209b. In detail, the second 
multiplication module 209a includes six second 
multiplication units 215 which multiply a respective second- 
layer output oI2 (n) , supplied by the second-layer output 
calculating unit 208, by a respective third-layer weight 
W DF {n) , supplied by the third work memory 203. The 
defuzzif ication module 209b calculates the reconstruction 
weights oL3, adding the products supplied by the second 
multiplication module 209a, adding together the second-layer 
outputs oL2 (n) , and calculating the ratio between the two 
siims in accordance with equation (6). 

The reconstruction unit 210 stores the reconstruction 
weights o!r3(i) as these are supplied by the reconstruction- 
weight calculating unit 209 and, as soon as it has stored 
2W + 1 reconstruction weights oL3, it calculates the output 
sample u{2N + 1) in accordance with equation (7), also using 
the values of the 2N + 1 samples supplied by. the data memory 
200. Subsequently, upon receipt of the next reconstruction 
weight oL3 (2N +2), it calculates the output sample 
u(2W + 2) using also the previous 2N reconstruction weights 
oL3 and as many samples coming from the data memory 200,. in 
accordance with equation (7) . 

The control unit 206 determines the processing sequence and 
data loading/transfer between the various units and modules. 
To calculate a single output sample, the control unit 
repeats the sequence of steps of the fuzzy algorithm 2N + 1 
times, updates the data memory 200 at the end of 2N + 1 
cycles, and controls loading of successive 2N + 1 samples. 

The advantages of the method and filter illustrated herein 
are described below. First, the method and filter reduce the 
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noise of the input signal, whether the noise in question is 
of a white type or of a colored type, and enable separation 
of signals having different features. The filter preserves 
the steep edges of the signals without causing any losses of 
5 signal features, as is evident from a comparison between 

Figures 10a, 10b, and 10c. In particular, Figure 10a is a 
plot of a non-noisy voice signal (namely, a signal fragment 
corresponding to the vowel "e", with sampling at 44.1 kHz 
and a 16-bit A-D conversion resolution) ; Figure 10b is a 
10 plot the same signal as in Figure 10a, in presence of white 

noise; and Figure 10c shows the result of filtration of the 
signal of Figure 10b using the filter device 95 of Figure 7. 

Since the filter can be trained, it can be adapted to a 
15 specific type of initial signal and can be subsequently 

modified if so required. For example, the filter can be 
initially adapted to a first type of acoustic signal (for 
instance, a male voice with white noise, training being 
performed with two signals, one in the presence of white 
20 noise and one in the absence of noise) . Then, the learning 

step can be repeated for signals of a different type, for 
example a female voice. The new weights thus calculated can 
replace the previous ones or be stored in addition to the 
previous ones; in the latter case, the control unit 206 
25 controlling the filter 1 can control sending, to the units 

207-209 (Figure 9), either of the first or of the second 
samples according to the use of the filter. 

Finally, it is clear that numerous variations and 
30 modifications may be made to the method and filter described 

and illustrated herein, all falling within the scope of the 
invention, as defined in the annexed claims. 
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- J CLAIMS 

1. A filtering device (80; 95) for filtering signals having 
steep edges, characterized by at least one neuro- fuzzy 
filter (1) . 

5 

2. A filtering device according to Claim 1, characterized in 
that said neuro-fuzzy filter (1) comprises: 

a signal-feature calculating unit (2) receiving input 
samples (e(i)) of a signal (In) to be filtered and 

10 generating signal features (XI (i), X2(i), X3(i)); 

a neuro-fuzzy network (3) receiving said signal features and 
generating reconstruction weights (oL3(i)); and 
a moving-average reconstruction unit (4) receiving said 
input samples (e(i)) and said reconstruction weights 

15 (oL3(i)), and generating output samples (u(i)) from said 
input samples and said reconstruction weights. 

3. A filtering device according to Claim 2, characterized in 
that said signal-feature calculating unit (2) comprises a 

20 memory (200) storing a first plurality of input samples 
(e(i)) forming a first window, and a feature-calculating 
network (201, 205) receiving said first plurality of input 
samples and supplying a plurality of features (XI (i) , X2(i), 
X3(i)) for each one of said input samples. 

25 

4. A device according to Claim 3, characterized in that said 
feature-calculating network (201, 205) comprises first 
feature-providing means (201) which generate, for each of 
said input samples (e(i)), a first signal feature (XI (i) ) 

30 correlated to a position of said input sample in said first 
window; second feature-providing' means (205a) which 
generate, for each of said input samples (e(i)), a second 
signal feature (X2(i)) correlated to the difference between 
said input sample and a central sample {e{N) ) in said first 

35 window; and third feature-providing means (205b) which 
generate, for each of said input samples (e(i)), a third 
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signal feature (X3 (i) ) correlated to the difference between 
said input sample and an average sample value in said first 
window . 

5 5. A device according to Claim 4, characterized in that said 
first feature-providing means (201) generate said first 
signal ' feature (XI (i) ) for an input sample (e(i)) according 
to the relation 
I/-JVI 

10 wherein i is the position of said input sample (e(i)) in 
said first window, and N is the position of a central sample 
in said first window. 

6. A device according to Claim 5, characterized in that said 
15 first feature-providing means comprise a first work memory 

(201) storing said first signal feature (XI (i) ) for each one 
of said input samples (e(i)). 

7. A device according to any of Claims 4-6, characterized in 
20 that said second feature-providing means (205a) generate 

said second signal feature (X2 (i) ) for an input sample 
(e(i)) according to the relation 

max(diff) 

wherein e(N) is a central sample in said first window, and 
25 max(di-f-f) is the maximum of the differences between all the 
input samples in said first window and said central sample. 

8. A device according to any of Claims 4-7, characterized in 
that said third feature-providing means (205b) generate said 

30 third signal feature (X3(i)) for an input sample (e(i)) 

according to the relation 

\e(j)-av\ 

X3(i)= 1 ' — 

max(diffav) 

wherein av is the average value of the input samples in said 
first window, and max ( diff_av) is the maximum of the 
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differences between all the input samples in said first 
window and said average value av. 

9. A device according to any of Claims 2-8, characterized in 
5 that said neuro-fuzzy network (3) comprises: 

fuzzif ication neurons (6, 202) receiving said signal 
features (XI (i), X2(i), X3 (i) ) of an input sample (e(i)) and 
generating first-layer outputs (oLl) defining a confidence 
level of said signal features with respect to preset 
10 membership functions; 

fuzzy neurons (7) of an AND type, receiving said first-layer 
outputs and generating second-layer outputs (oL2) deriving 
from fuzzy rules; and 

a defuzzif ication neuron (8) receiving said second-layer 
15 outputs and generating a reconstruction weight (oL3) for 
each of said input samples (e(i)), using a center-of-gravity 
criterion . 

10. A device according to Claim 9, characterized in that 
2 0 said membership functions are gaussian functions, and said 

first-layer outputs are calculated according to the equation 

' (Xl-W m (l,k)~ 



oLl(l,k) = exp 



wherein oil {I, k) is a first-layer output, XI is a signal 
feature, W m {l, k) is the mean value, and W V {1, k) is the 
25 variance of a gaussian function. 

11. A filtering device according to Claim 9 or 10, 
characterized in that it comprises two membership functions 
for each one of said signal features (XI (i) , X2(i), X3(i)). 

30 

12. A filtering device according to any of Claims 9-11, 
characterized in that said first-layer neurons (6) comprise 
a second work memory (202) storing values of said first- 
layer outputs (oil (1, k) ) for each value of said signal 

35 features (XI (i) , X2(i), X3(i)). 
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13. A filtering device according to any of Claims 9-12, 
characterized in that said fuzzy rules are of an AND type, 
and said second-layer outputs (oL2) use the norm of the 

5 minimum. 

14. A filtering device according to Claim 13, characterized 
in that said second-layer outputs are calculated according 

to the equation 
10 oL2(n) = mm§V FA (m,n)-oLl(m)} 

n 

wherein oL2 (n) is a second-layer output; W rA (m, n) is a 
second-layer weight, and oLl (1, k) is a first-layer output. 

15. A filtering device according to Claim 14, characterized 
15 in that said second-layer neurons (7) comprise a plurality 

of multiplication units (213) and a plurality of minimum 
modules (214) cascade-connected together. 



16. A filtering device according to any of Claims 9-15, 
characterized in that said reconstruction weights are 
calculated according to the equation 
^W DF {ri)-oL2{ri) 



oL3=— 



wherein oL3 (n) is a reconstruction weight; W DF {n) are third- 
layer weights, and oL2 (n) is a second-layer output. 



17. A filtering device according to any of Claims 9-16, 
characterized in that said moving-average reconstruction 
unit (4) receives a second plurality of input samples (e(i)) 
forming a second window, and a corresponding plurality of 
30 reconstruction weights (o£3(i)), and calculates each of said 
output samples (u(i)) according to the equation 



4 
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2N 

^oL3(i-j).e(i-j) 

«(0= i=1 -5v 

2>('-J) 

1=0 

wherein e(i-j) is an (i— j)-th input sample, and oI/3(i-j) is 
a reconstruction weight associated to an (i— j)-th input 
sample ." 

18. A filtering device according to any of Claims 9-17, 
characterized by a training unit (5) having a first input 
connected to said moving- average reconstruction unit (4) and 
receiving said output samples (u(i)), a second input 
receiving a desired output signal ( Tg) , and an output 
connected to said neuro-fuzzy network (3) to supply 
optimized weighting values. 

19. A filtering device according to any of Claims 1-18, 
characterized by a first splitting stage [Hi, H 0 , 84, 85; 
•Hooi, Hi 02, 96, 97) generating at least two streams of samples 
to be filtered; one said neuro-fuzzy filter (86, 87; 120- 
127) for each stream of samples to be filtered, each of said 
neuro-fuzzy filters generating a respective stream of 
filtered samples; and a first recomposition stage (88, 89, 
92, G 0 , Gi; 154, 155, 160, Goon G102) receiving said streams 
of filtered samples and generating a single stream of output 
samples (s 0 ) • 

20. A filtering device according to Claim 19, characterized 
by a plurality of further splitting stages (H011-H128/ 100- 
114) cascade-connected together and to said first splitting 
stage (H001, #102, 96, 97), and a plurality of further 
recomposition stages (130-151, G021-G114) cascade-connected to 
each other between said neuro-fuzzy filters (120-127) and 
said first recomposition stage (154, 155, 160, G 00 i, G102) . 

21. A filtering device according to Claim 19 or 20, 
characterized in that said splitting stages [Hi, H 0 , 84, 85; 



Hooi-#i23 / 96-114) each comprise a first and a second analysis 
filters (Hi, Ho, Hooi-Hiza) in phase quadrature to each other 
and receiving a stream of samples to be split, said first 
and a second analysis filters generating a respected stream 
of split samples, and a first and a second downsampling unit 
(84, 85; 96-114), each of which receives a respective stream 
of split samples, 

and in that said recomposition stages (88, 89, 92, G 0 , G x ; 
130-160, G021-G102) each comprise a first and a second 
upsampling units (88, 89; 130-137, 146-149, 154, 155), each 
first and a second upsampling units receiving a respective 
stream of samples to be incremented and generating a 
respective stream of incremented samples; a first and a 
second synthesis filters {Go, d; G 0 :i-Gi 02 ) in quadrature with 
to each other and complementary to said analysis filters, 
each of said first and a second synthesis filters receiving 
a respective stream of incremented samples and generating a 
respective stream of partial samples; and an adder node (92; 
140-143, 150, 151, 160) receiving said streams of partial 
samples and generating a stream of added samples. 

22. A filtering device according to Claim 21, characterized 
in that said analysis filters (Hi, H 0 ; H001-H129) are 
quadrature mirror filters (QMFs) , and said synthesis filters 
(G 0 , Gi; G021-G102) are QMFs complementary to said analysis 
filters . 

23. A filtering device according to Claim 22, characterized 
in that said quadrature mirror filters (Hi, H 0 ; H 0 oi-Hi2 8 ; G 0 , 
Gi; G021-G102) are convolutive filters. 

24. A method for reducing noise in a filter having sharp 
edges, characterized by a neuro-fuzzy filtering step. 



25. A method according to Claim 24, characterized in that it 
comprises the steps of: 



calculating signal features (XI (i) , X2(i), X3(i)) from input 
samples (e(i)) of a signal {In) to be filtered; 

calculating reconstruction weights (oL3) from said signal 
features using a neuro-fuzzy network (3) ; and 

reconstructing, from said input samples and said 
reconstruction weights and using a moving-average filter, an 
output signal (Out) including a plurality of output samples 
(uU) ) • 

2 6. A method according to Claim 25, characterized by the 
steps of: 

storing a first plurality of input samples (e(i)) forming a 
first window; and 

calculating, from said first plurality of input samples, a 
plurality of signal features (XI (i), X2(i), X3(i)) for each 
of said input samples. 

27. A method according to Claim 26, characterized in that 
said step of calculating a plurality of signal features for 
each of said input samples comprises the steps of: 
calculating a first signal feature (XI (i)) correlated to a 
position of said input sample (e(i)) in said first window; 
calculating a second signal feature (X2 (i) ) correlated to 
the difference between said input sample and a central 
sample in said first window; and 

calculating a third signal feature (X3 (i) ) correlated to the 
difference between said input sample and an average sample 
value av in said first window. 

28. A method according to Claim 26 or 27, characterized in 
that said step of calculating reconstruction weights (oL3) 
comprises the steps of: 

performing a fuzzif ication operation and calculating first- 
layer outputs (oil) defining confidence levels of said 
signal features (XI (i), X2(i), X3(i)) with respect to preset 
membership functions; 



oLl(/,fc)=exp 



V 



performing a fuzzy AND operation and generating second-layer 
outputs (oL2) deriving from fuzzy rules, starting from said 
first-layer outputs; and 

performing a defuzzif ication operation on said second-layer 
outputs and generating a reconstruction weight (o£3) for 
each one of said input samples (e(i)), using a center-of- 
gravity criterion of the. 

29. A method according to Claim 28, characterized in that 
said membership functions are gaussian functions, and said 
first-layer outputs are calculated according to the equation 

' Xl-W m (l,k) 

WM,*) . 

wherein oil (1, k) is a first-layer output, XI is a signal 
feature, W m (l, k) is the mean value, and W v (l, k) is the 
variance of a gaussian function. 

30. A method according to Claim 28 or 29, characterized in 
that said second-layer outputs are calculated according to 

the equation 

oL2(ri) = min §V FA (m, w) • oL\(m)} 

n 

wherein oL2(n) is a second-layer output; W FA (m, n) is a 
second-layer weight, and oLIU, k) is a first-layer output. 

31. A method according to any of Claims 28-30, characterized 
in that said reconstruction weights are calculated according 
to the equation 

jrW DF (ri)oL2(ri) 



oL3=~ 



2>£2(«) 

n=l 

wherein oL3 is a reconstruction weight; W DF (n) are third- 
layer weights, and oL2 (n) is a second-layer output. 

32. A method according to any of Claims 26-31, characterized 
in that said step of reconstructing using a moving-average 
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filter comprises the steps of: 

receiving a second plurality of input samples (e(i)) forming 
a second window, and a respective plurality of 
reconstruction weights (oL3(i)); and 

calculating each of said output samples according to the 
equation 

2N 

^oL3{i-j)-e(i-j) 

^>^—r N 

wherein e(i-j') is an (i-j')-th input sample, and oL3(i—j) is 
a reconstruction weight associated to an (i— j')-th input 
s amp 1 e . 

33. A method according to any of Claims 24-32, characterized 
by a training step for training weights used in said neuro- 
fuzzy filtering step. 

34. A method according to Claim 33, characterized in that 
said training step comprises the steps of: 

generating an input signal having a known configuration; 
filtering said input signal having a known configuration to 
obtain a test output signal; 

comparing said test output signal with a desired signal ( Tg) 
to obtain a distance between said test output signal and 
said desired signal; 

calculating a fitness function from said distance; and 
optimizing said weights in accordance with said fitness 
function . 

35. A method according to any of Claims 2 6-34, characterized 
by a multiresolution analysis whereby the signal is splitted 
into sub-bands through orthonormal wavelets . 

36. A method according to Claim 35, characterized by the 
steps of: 

splitting a stream of input samples into at least two 
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streams of samples to be filtered; 

filtering each stream of samples to be filtered using a 
respective neuro-fuzzy filter (86, 87; 120-127) to obtain at 
least two streams of filtered samples; and 
5 recomposing said streams of filtered samples to generate a 
single stream of output samples (s c ) . 

37. A method according to Claim 36, characterized in that, 
before performing said step of filtering, said step of 

10 splitting is repeated a preset number of times, and in that, 
after performing said step of filtering, said step of 
recomposing is repeated said preset number of times. 

38. A method according to Claim 36 or 37, characterized in 
15 that said step of splitting comprises the steps of: 

feeding a stream of samples to be split to two analysis 
filters in phase guadrature to each other; 
generating two streams of filtered split samples; and 
downsampling said streams of filtered split samples, 
20 and in that said step of recomposing comprises the steps of: 
upsampling streams of samples to be incremented, generating 
streams of incremented samples; 

filtering said streams of incremented samples using two 
synthesis filters in phase quadrature to each other and 
25 complementary to said analysis filters, generating streams 
of partial samples; and 

adding pairs of streams of partial samples and generating a 
stream of added samples. 
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ABSTRACT 



The filtering device (80) comprises a neuro-fuzzy filter (1; 
80) and implements a moving- average filtering technique in 
which the weights for final reconstruction of the signal 
5 (oI3(i)) are calculated in a neuro-fuzzy network (3) 
according to specific fuzzy rules. The fuzzy rules operate 
on three signal features (XI (i) , X2 (i) , X3(i)) for each 
input sample (e(i)). The signal features are correlated to 
the position of the sample in the considered sample window, 

10 to the difference between a sample and the sample at the 
center of the window, and to the difference between a sample 
and the average of the samples in the window. The filter 
device for the analysis of a voice signal comprises a bank 
of neuro-fuzzy filters (86, 87) . The signal is split into a 

15 number of sub-bands, according to wavelet theory, using a 
bank of analysis filters including a pair of FIR QMFs {H 0 , 
Hi) and a pair of downsamplers (85, 86); each sub-band signal 
is filtered by a neuro-fuzzy filter (86, 87), and then the 
various sub-bands are reconstructed by a bank of synthesis 

20 filters including a pair of upsamplers (88, 89), a pair of 
FIR QMFs (Go, Gi), and an adder node (92). 

Figures 1 and 4 
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IZflMAGE CUT OFF AT TOP, BOTTOM OR SIDES 
fcj FADED TEXT OR DRAWING 
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□ GRAY SCALE DOCUMENTS 
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